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ABSTRACT 


This  report  describes  the  EVD  Jet-Wing  Computer  Program,  which  is  based 
upon  the  Elementary  Vortex  Distribution  (EVD)  Jet-Wing  Lifting  Surface  Theory 
described  in  Volume  I of  this  report.  This  program  provides  a capability  for 
determining  the  aerodynamic  characteristics  of  wings  of  arbitrary  planform, 
and  includes  the  following: 

1.  Spanwise  and  chordwise  loading 

2.  Spanwise  variation  of  induced  drag 

3.  A capability  to  investigate  the  effects  of: 

a.  Part  span  flaps 

b.  Part  span  blowing 

c.  Pitching,  rolling,  yawing  and  sideslip 

4.  Total  lift  and  induced  drag  (momentum  method),  pitching,  yawing 

and  rolling  moments,  etc. 

The  program  has  the  capabilities  for  investigating  the  effects  of  a 
variation  of  leading  and  trailing  flap  deflection,  camber,  twist,  jet  deflec- 
tion and  jet  momentum. 
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PREFACE  TO  REVISED  EDITION 


This  report  was  originally  published  in  May  1973  in  two  volumes:  Volume  I, 
The  Elementary  Vortex  Distribution  Jet-Wing  Lifting  Surface  Theory  and 
Volume  II,  EVD  Jet-Wing  Computer  Program  User's  Manual.  In  this  new  edition 
several  typographical  errors  have  been  corrected,  a minor  change  made  In  the 
computer  program,  and  the  "Limited  Distribution"  of  Volume  II  removed. 


Michael  L.  Lopez 
Douglas  Aircraft  Company 
McDonnell  Douglas  Corporation 
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NOMENCLATURE 


INPUT/ 


NAME 

SYMBOL 

OUTPUT 

DEFINITION 

A 

A 

I/U 

Fundamental  case  scale  factor 

AC 

Cc 

1/0 

Local  incidence  due  to  camber 

ACTE 

ect 

1/0 

Local  incidence  due  to  camber  at  trailing  edge 

ALFIN 

0 

Downwash  angle  or  jet  angle  at  infinity 
downstream 

ALPHA 

a 

0 

Angle  of  attack 

ARATI0 

AR 

1/0 

Aspect  ratio 

AREA 

S 

1/0 

Wing  planform  reference  area 

BETA 

B 

1/0 

Deflection  angle,  (6^,  6j,  or  6$) 

CBGL 

V 

0 

Root  bending  moment  coefficient  due  to 
pressure  on  left  wing 

CBGR 

°BrR 

0 

Root  bending  moment  coefficient  due  to 
pressure  on  right  wing 

CBUL 

c8jl 

0 

Root  bending  moment  coefficient  due  to 
jet  reaction  on  left  wing 

CBJR 

C“Jr 

0 

Root  bending  moment  coefficient  due  to 
jet  reaction  on  right  wing 

CBL 

cbl 

0 

Total  root  bending  moment  coefficient 
on  left  wing 

CBR 

cBr 

0 

Total  root  bending  moment  coefficient 
on  right  wing 

CCD 

(Cd1}P 

0 

Total  induced  drag  coefficient  calculated 
by  pressure  integration 

CCJ 

CJ 

0 

Total  jet  momentum  coefficient 

CCL 

CL 

0 

Total  lift  coefficient 

CCM 

^m 

0 

Total  pitching  moment  coefficient  (about 
wing  apex) 

CCS 

CS 

0 

Total  leading  edge  suction  coefficient 

lx 


NOMENCLATURE 


INPUT/ 


NAME 

SYMBOL 

OUTPUT 

DEFINITION 

CCT 

CT 

0 

Total  net  thrust  coefficient 

CCY 

CY 

0 

Total  side  force  coefficient 

CD 

cdi 

0 

Sectional  induced  drag  coefficient 

CDITZ 

(cDi)M 

0 

Total  induced  drag  coefficient,  calculated 
by  the  momentum  method 

CH0RD 

c 

0 

Wing  sectional  chord 

CL 

c* 

0 

Sectional  lift  coefficient 

CLL 

Ci 

0 

Rolling  moment  coefficient 

CLLP 

% 

0 

Rolling  moment  coefficient  derivative  due 
to  rolling 

CLLR 

C'r 

0 

Rolling  moment  coefficient  derivative  due 
to  yawing 

CLQ 

% 

0 

Lift  coefficient  derivative  due  to  pitching 

CM 

cm 

0 

Sectional  pitching  moment  coefficient 
(about  the  local  leading  edge) 

CMAC 

c 

0 

Mean  aerodynamic  chord  (MAC) 

CMQ 

Cmq 

0 

Pitching  moment  coefficient  derivative  due 
to  pitching  about  center  of  gravity 

CMU 

c 

u 

1/0 

Sectional  jet  momentum  coefficient 

CN 

Cn 

0 

Yawing  moment  coefficient 

CNI 

C"i 

0 

Total  yawing  moment  coefficient  due  to 
induced  drag  (including  leading  edge  suction) 

CN(P) 

CnP 

0 

Yawing  moment  coefficient  derivative  due 
to  rolling,  dependent  on  rolling  rate 

CN(R) 

cnr 

0 

Yawing  moment  coefficient  derivative  due 
to  yawing,  dependent  on  yawing  rate 

CP 

*CP 

0 

Pressure  coefficient,  Acp=cpi()Wer  - cPupper 

surface  surface 

x 
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NOMENCLATURE 

NAME 

SYMBOL 

INPUT/ 

OUTPUT 

DEFINITION 

CPMBL 

yc.l. 

wr- 

0 

Spanwise  center  of  lift  of  left  wing,  in 
terms  of  half  span 

CPMBR 

Yc.l. 

W 

0 

Spanwise  center  of  lift  of  right  wing,  in 
terms  of  half  span 

CREF 

c 

1/0 

Wing  reference  chord 

CS 

cs 

0 

Sectional  leading  edge  suction  coefficient 

CTO 

ct 

0 

1.  Sectional  thrust  coefficient 

2.  Tip  chord  of  wing 

CXCL 

CXCLB 

XC.L. 

C 

XC.L. 

B72“ 

0 

0 

x-coordinate  of  center  of  lift,  in  terms 
of  reference  chord 

CXCL  = - (CCMG  + CCML)/CCL  labeled  (X/CREF) 

x-coordinate  of  center  of  lift,  in  terms 
of  half  span 

CXCLB  = - ((CCMG  + CCML)/CCLJ* CREF  labeled 
(X/B/2) 

CXCP 

XC.P. 

c 

0 

x-coordinate  of  center  of  pressure,  in 
te^ns  of  reference  chord 

CXCPB 

XC.P. 

~B7 r 

0 

x-coordinate  of  center  of  pressure,  in 
terms  of  half  span 

CY(P) 

% 

0 

Side  force  coefficient  derivative  due  to 
rolling,  dependent  on  rolling  rate 

CY  ( R) 

CYr 

0 

Side  force  coefficient  derivative  due  to 
yawing,  dependent  on  yawing  rate 

D 

d 

0 

Chordwise  distance  from  trailing  edge  to 
infinity  EVD  vortex  point 

DEL 

? 

0 

Chordwise  length  of  an  EVD  element  in 
terms  of  local  chord 

DELTA 

A 

0 

Half  of  the  spanwise  width  of  an  EVD 
element  or  a spanwise  division,  normalized 
by  b/2 

DJ 

6J 

1/0 

Jet  deflection  angle  relative  to  the  wing 
trailing  edge  slope 

NOMENCLATURE 


INPUT/ 


NAME 

SYMBOL 

OUTPUT 

DEFINITION 

EPS 

ei 

0 

Total  local  incidence  angle  at  a downwash 
control  point  on  the  wing 

GAMMA 

Y 

0 

Vorticity  intensity 

HL 

\ 

1/0 

Height  of  leading  edge  above  the  x-y  plane 

I 

i 

0 

EVD  element  sequence  number 

ICT 

I 

Camber  type  flag  for  each  wing  section 

ICTYPE 

I 

Wing  chordwise  division  type  flag  for  each 
section 

IDERIV 

I 

Stability  derivative  flag 

IGTYPE 

1/0 

Planform  geometry  type  flag 

IHINGE 

1/0 

Hinge  EVD  flag 

IHT 

I 

Hinge  type  flag  for  each  wing  section 

IJTYPE 

I 

Jet  chordwise  division  type  flag  for  each 
section 

ILT 

I 

Flap  type  flag 

INBETA 

I 

Hinge  angle  input  flag 

INCAMB 

I 

Camber  angle  input  flag 

INDELJ 

I 

Jet  angle  input  flag 

INHITE 

I 

Leading  edge  height  input  flag 

INTWST 

I 

Twist  angle  input  flag 

IPRINT 

1/0 

Output  control  flag 

ISYMM 

1/0 

Symmetry  control  flag 

J 

J 

Sectional  jet  momentum 

JETFLG 

1/0 

Jet  control  flag 

K 

k 

Sequence  number  of  spanwise  divisions  or 
sections 

NOMENCLATURE 


I 


NAME 

SYMBOL 

INPUT/ 

OUTPUT 

N 

1/0 

NCASES 

K 

1/0 

NCT 

I 

NHT 

I 

NI 

I 

NJ 

0 

NJTYPE 

I 

NR0WS 

0 

NR0WSJ 

I 

NW 

0 

NWTYPE 

I 

SPAN 

b 

1/0 

SWEEP 

Ac/4 

1/0 

TANLE 

tanA^ 

0 

THETA 

e 

0 

TITLE 

1/0 

TR 

A 

I 

TWIST 

“o 

1/0 

TYPE 

0 

U 

U 

X 

X 

0 

DEFINITION 

Fundamental  case  sequence  number 

Number  of  fundamental  rases 

Number  of  camber  types 

Number  of  hinge  types 

Number  of  chordwise  divisions  of  a 
particular  chordwise  division  type  (wing 
or  jet) 

Number  of  chordwise  divisions  of  jet  section 
Number  of  jet  chordwise  division  types 
Number  of  wing  sections  or  spanwise  divisions 
Number  of  sections  having  a trailing  jet 
Number  of  chordwise  division  of  each  section 
Number  of  wing  chordwise  division  types 
Wing  span 

Wing  quarter-chord  sv/eep  angle 

Tangent  of  sectional  leading  edge  sweep  angle 

Total  jet  deflection  angle  relative  to  the 
free  stream 

Run  title 

Wing  taper  ratio 

Twist  angle  for  each  wing  section 

EVD  type  flag 

Free  stream  velocity 

Chordwise  coordinate  of  the  Cartesian 
coordinate  system 
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I 


xiii 


NAME 

XB 

XBH 

XCL/C 

XCP/C 

XCG 

XI 

XLEAD 

XMC 

XTRAIL 

Y 


SUFFIX 

A 

A2 

G 

J 

MC 

Mil 


NOMENCLATURE  ~ — 


INPUT/ 
SYMBOL  OUTPUT 

X 1/9 

xh  1/9 


c.g. 


1/9 


xi  9 

x,  1/(3 

Vc.  w 


H W 

y 1/(3 


DEFINITION 

Chordwlse  distance  of  an  EVD  element  from 
the  leading  edge,  In  terms  of  local  chord 

Chordwise  distance  of  a hinge  point  from 
the  leading  edge,  In  terms  of  local  chord 

Chordwlse  location  of  center  of  lift  at  a 
section  measured  from  local  leading  edge, 
in  terms  of  local  chord 

Chordwise  location  of  center  of  pressure 
at  a section  measured  from  local  leading 
edge,  in  terms  of  local  chord 

x-coordinate  of  the  center  of  gravity  about 
which  pitching  rate  derivatives  are  taken 

x-coordinate  of  an  EVD  element 

x-coordinate  of  the  leading  edge 

x-coordinate  of  the  moment  center  about 
which  pitching  moments  are  taken 

x-coordinate  of  the  trailing  edge 

y-coord.inate  of  a section 


Subscripts 


SUBSCRIPT  DEFINITION 

a Indicates  a linear  variation  with  angle  of  attack 

2 

<*  Indicates  a quadratic  variation  with  angle  of  attack 

r Indicates  a contribution  due  to  pressure  (circulation) 

J Indicates  a contribution  due  to  jet  reaction 

m.c.  Indicates  pitching  moments  taken  about  moment  center 

u Indicates  a contribution  due  to  jet  reaction  at  a wing 

section 
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NOMENCLATURE 


SUFFIX 

SUBSCRIPT 

DEFINITION 

R2 

r2 

Indicates  a quadratic  variation  with  yawing  rate 

T 

t 

Indicates  a contribution  due  to  the  thrust  component 
of  the  jet  reaction  at  a section 

X 

a 

Indicates  a contribution  from  mutual  Interference 
between  the  basic  configuration  and  the  angle  of 

attack 

0 

0 

Indicates  a * 0 
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1.0  INTRODUCTION 


This  report  describes  the  Mark  II  version  of  the  EVD  Jet-Wing  Computer 
Program,  intended  for  use  as  a preliminary  design  tool  for  the  analysis  of 
arbitrary  wings  incorporating  trailing  edge  blowing.  The  program  is  based 
upon  the  linearized  Elementary  Vortex  Distribution  (EVD)  Jet-Wing  Lifting 
Surface  Theory  described  in  Volume  I of  this  report  (Reference  1).  The 
program  is  written  in  the  Fortran  IV  language  and  can  be  adapted  for  use  on 
many  existing  large  computing  sy terns.  Versions  are  currently  in  operation  on 
the  CDC  6000  series  and  IBM  360  and  370  series  computers. 

In  order  to  execute  the  program,  the  user  must  describe  the  arbitrary 
geometric  configuration  to  the  computer  by  using  several  different  types  of 
input  cards  containing  specific  information  in  a fixed  order.  Deflection,  for 
example,  of  a wing  control  surface  may  be  accomplished  by  using  the  linear 
"fundamental  case"  feature  of  the  program.  Superposition  of  the  aerodynamic 
characteristics  associated  with  fundamental  cases  so  as  to  provide  the  overall 
aerodynamic  characteristics  of  a complete  configuration  may  be  obtained  by 
specifying  the  "composite  case"  feature  of  the  program.  The  number  of  jet 
strength  cases  which  can  be  produced  for  each  geometric  configuration  is 
limited  only  by  the  computing  time  available.  Printed  output  from  the  program 
contains  the  jet-wing  aerodynamic  characteristics  in  coefficient  form.  Detailed 
spanwise  and  chordwise  loading  and  certain  dynamic  stability  derivatives  may 
also  be  included  at  the  option  of  the  user. 

Maximum  realization  of  the  potential  of  this  program  can  only  be  achieved 
by  a user  who  is  thoroughly  familiar  with  the  applications  and  limitations  of 
the  linearized  theory  upon  which  it  is  based,  and  who  is  willing  to  apply  his 
own  engineering  knowledge  and  judgment  in  a methodical  manner  to  each  con- 
figuration analyzed.  It  is  suggested  that  both  volumes  of  this  report  be 
read  before  running  the  first  problem,  and  that  a few  simple  problems  be 
tried  in  order  to  develop  understanding  and  familiarity  with  the  program 
before  large  complicated  problems  are  attempted. 


2.0  FUNDAMENTALS  OF  PROGRAM  USE 


The  EVD  jet-wing  computer  program  is  based  upon  the  linearized  EVD 
jet-wing  lifting  surface  theory  described  in  Reference  1.  This  program  has 
several  unique  features  which  enable  it  to  be  used  by  the  aerodynamicist  in 
the  evaluation,  design  and  development  of  V/STOL  aircraft  systems  which  utilize 
powered  high  lift  systems  based  on  the  jet  flap  principle  (e.g.,  EBF,  augmentor 
wing,  etc.).  The  scope  of  this  program  is  not  limited,  however,  to  the  analysis 
of  jet  flap  wings.  By  selecting  several  options  available  to  the  user,  this 
program  may  be  used  as  a classical  lifting  surface  method.  Hence,  as  an  aero- 
dynamic tool,  the  analysis  of  conventional  wings  can  be  regarded  as  being 
within  the  scope  of  the  present  computer  program. 

The  program  is  written  entirely  in  the  Fortran  IV  language,  and  can  be 
adopted  for  use  on  any  large-scale  computer  system.  The  program  consists  of 
several  primary  components  designed  to  minimize  computer  storage  requirements 
while  maintaining  maximum  computer  efficiency  and  speed.  This  modular  arrange- 
ment also  results  in  easier  program  checkout  and  correction. 

The  programming  philosophy  has  sought  to  maintain  the  independence  of 
each  small  component  of  the  program  wherever  possible.  This  approach  also 
facilitates  the  addition  of  new  capabilities,  with  minimum  alteration  of  the 
existing  components. 

2.1  Philosophy  of  Program  Use 

It  is  fundamentally  important  that  the  user  of  the  EVD  jet-wing  computer 
program  understand  the  basic  limitations  and  approximations  adopted  in  the 
solution  of  the  jet-vn'ng  problem.  The  program  is  not  a "magic  box"  with  which 
he  can  expect  to  achieve  valid  results  while  supplying  poorly  prepared  input 
data.  He  must,  therefore,  understand  each  step  of  his  role  in  the  analysis 
of  a specific  problem,  including  use  of  engineering  knowledge  and  judgement  in 
every  decision  he  makes  in  the  preparation  of  computer  input  data.  Only  then 
can  he  expect  the  program,  which  is  essentially  an  engineering  tool,  to  reliably 
and  accurately  complete  the  task  assigned  to  it.  This  point  may  be  summarize 
by  the  quaint  term  which  applies  to  all  computer  programs:  "GIGO"  - Garbage 
In,  Garbage  OutJ 
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Judgement  is  also  required  in  the  interpretation  and  application  of  the 
computer  output  data.  The  user  should  be  aware  of  those  factors  not  included 
in  the  program  analysis,  which  may  have  an  effect  on  the  final  result.  He 
should  also  have  some  understanding  of  the  way  in  which  these  factors  would  be 
likely  to  influence  the  program  results.  For  example,  if  there  are  regions  of 
flow  separation  on  a real  wing,  the  program  would  be  expected  to  overpredict 
the  lift  of  that  wing.  On  the  other  hand,  a thick  wing  is  expected  to  have 
higher  lift  than  the  wing  of  zero  thickness  which  is  considered  by  the  present 
program.  The  quantitative  effects  of  such  factors  are,  of  course,  very  diffi- 
cult to  accurately  predict,  but  the  user  should  at  least  be  aware  of  the  trends 
which  they  are  expected  to  produce. 

The  user  is  also  cautioned  against  extrapolating  the  linearized  results 
into  regions  where  strong  nonlinear  behavior  might  be  expected.  For  example, 
the  small  angle  assumptions  of  the  linearized  approach  make  it  unlikely  that 
the  program  would  accurately  predict  the  characteristics  of  a wing  with,  say, 
a flap  deflected  at  60  degrees.  While  the  program  would  give  an  indication  of 
the  trends  produced  by  such  a nonlinear  configuration,  and  in  certain  cases 
linearized  results  have  been  shown  to  be  remarkably  good  in  seemingly  nonlinear 
situations,  it  would  nevertheless  be  dangerous  to  use  the  output  data  as  though 
it  were  generally  an  accurate  prediction  of  the  aerodynamic  characteristics 
under  all  conditions. 

2.2  Brief  Review  of  the  Linearized  Approach 

The  assumptions  and  restrictions  of  the  EVD  jet-wing  lifting  surface 
theory  are  discussed  in  detail  in  Volume  I of  this  report.  For  the  convenience 
of  the  reader,  the  basic  assumptions  are  restated  here: 

On  the  Wing: 

a.  The  wing  is  thin,  and  is  represented  by  the  mean  camber  line. 

b.  All  local  incidences  are  small,  but  may  be  discontinuous. 

On  the  Jet: 

a.  The  jet  is  thin. 

b.  The  jet  deflection,  relative  to  the  freestream,  is  small. 

General  Flowfield: 

a.  Incompressible,  Inviscid,  irrotational . 

b.  No  mixing  occurs  between  the  jet  and  the  external  flow. 


c.  All  spanwise  components  of  velocity  are  considered  small 
relative  to  the  freestream  velocity. 

d.  Rollup  of  the  jet  sheet  and  wing  and  jet  trailing 
vortex  system  is  neglected. 

The  above  assumptions  imply  the  existence  of  certain  other  conditions 
which  are  of  practical  importance  in  the  use  of  the  program.  Though  they  may 
not  seem  to  be  obviously  derived  from  the  basic  assumptions,  these  conditions 
nevertheless  are  controlling  factors  in  the  application  of  the  method,  and 
must  be  recognized  and  considered  by  the  user.  Such  conditions  include  the 
following: 

a.  The  wing  and  jet  are  represented  by  planar  sheets  of  bound 
trailing  vorticity. 

b.  All  potential  field  influences  of  the  wing  and  jet  originate 
in  the  wing  and  jet  plane. 

c.  All  boundary  conditions  are  satisfied  in  the  wing  and  jet 
plane. 

d.  The  jet  issues  from  the  trailing  edge  of  the  wing,  and 
its  reaction  force  acts  there  also. 

e.  Since  the  lifting  surface  theory  is  for  a jet-wing  system 
alone,  no  account  is  taken  for  the  effects  of  fuselage, 
nacelles,  empennage,  etc. 

f.  There  are  no  losses  or  dissipation  associated  with  the  jet 
sheet  emergence  from  the  wing  trailing  edge. 

g.  There  are  no  gaps  or  slots  in  the  wing  (program  restriction 
only). 

2 . 3 Geometry  Preparation 

The  EVD  jet-wing  computer  program  is  capable  of  evaluating  the  aerody- 
namic characteristics  of  jet-wings  of  arbitrary  wing  planform  and  arbitrary 
jet  spanwise  distribution  of  momentum.  In  all  cases,  the  jet  sheet  correctly 
extends  to  infinity  downstream.  The  jet-wing  system  may  be  considered  as 
symmetric,  anti -symmetric,  or  non-symmetric.  A symnetric  jet-wing  is  defined 
as  one  for  which,  under  all  conditions,  the  computed  aerodynamic  loading  would 
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be  symmetric  about  the  x-axis.  This  implies,  therefore,  symmetry  about  the 
x-axis  of  all  the  following  characteristics: 

Wing: 

a.  Planform  geometry 

b.  Local  incidence  of  each  element 

Jet: 

c.  Spanwise  location 

d.  Strength  at  each  spanwise  section 

e.  Deflection 

Anti-symnetric  jet-wings  are  defined  as  having  anti -symmetric  computed 
aerodynamic  loading.  Thus,  they  must  be  symmetric  in  all  of  the  above  items 
except  for  (b)  and  (e),  which  must  be  anti-symmetric.  Since  anti-symmetry  in 
the  angle  of  attack  of  the  wing  is  meaningless,  the  program  will  always  treat 
the  anti -symmetric  jet-wing  as  being  at  zero  angle  of  attack  only.  If  an 
otherwise  anti -symmetric  wing  were  at  an  angle  of  attack,  its  loading  would 
be  non- symmetric,  and  it  would,  therefore,  have  to  be  treated  as  completely 
non-symmetric.  Non-symmetric  jet-wings  are  defined  as  having  non-symmetry 
in  any  one  or  more  of  the  characteristics  listed  in  Items  (a)  through  (e) 
above. 


It  is,  however,  possible  for  the  user  to  superimpose  by  hand  the  results 
of  anti-symnetric  runs  with  those  of  symmetric  or  non-symmetric  runs.  Such  a 
situation  may  arise,  for  example,  when  a quantity  of  data  has  been  generated 
for  various  configurations  of  a particular  wing  and  the  effects  of  anti- 
symmetric ailerons  are  then  required.  It  may  be  more  efficient  to  compute  the 
aileron  effects  alone  in  one  additional  run  rather  than  rerunning  all  the 
previous  configurations  non-symmetrically  with  ailerons  included.  The  super- 
imposed loading  will  be  in  general  non-symmetric  and  both  the  right  and  left 
halves  of  the  wing  must  be  considered.  The  linear  aerodynamic  coefficients 
(e.g.,  CL,  CM,  Cj)  may  be  superimposed  directly  at  zero  angle  of  attack,  and 
the  variation  with  angle  of  attack  will  be  the  original  symmetric  or  non- 
symmetric  values.  The  non-linear  coefficients  (e.g.,  Cp,  Cy,  Cn)  cannot  be 
superimposed  directly.  Instead,  they  must  be  obtained  by  an  integration  of 
the  chordwise  and  spanwise  loading. 
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The  problem  analysis  begins  with  the  user  definition  of  the  planform 
geometry  of  the  wing  and  jet.  The  wing  and  jet  are  broken  up  into  spanwise 
sections  of  arbitrary  width.  The  arrangement  of  these  sections  and  their 
size  are  mainly  determined  by  two  factors.  First  is  the  occurrence  of  plan- 
form  discontinuities;  for  example,  flap  edges,  flap  extensions,  leading  and 
trailing  edge  brakes,  jet  limits,  etc.  Second  is  the  allowance  for  spanwise 
variation  of  loading.  Since  the  program  loading  is  assumed  to  be  that  which 
occurs  at  the  center  of  each  section  and  is  constant  in  the  spanwise  direction 
on  each  section,  there  should  be  close  spanwise  spacing  in  regions  where  rapid 
loading  variation  might  he  expected.  For  example,  near  the  wing  tip  or  near 
the  edge  of  a highly  loaded  flap  or  strong  jet,  there  will  normally  be  a rapid 
change  in  spanwise  loading.  The  spanwise  sections  should,  therefore,  be  narrow 
in  these  regions  so  that  the  "stairstep"  loading  of  the  program  can  adequately 
represent  a smooth  loading  variation. 

Each  spanwise  section  is  further  broken  up  into  chordwise  divisions  so 
that  the  jet-wing  planform  is  represented  by  an  array  of  rectangular  elements 
with  which  are  associated  an  equal  number  of  Elementary  Vortex  Distribution 
elements.  This  representation  is  shown  in  Figure  1 for  a typical  jet-wing. 

Note  that  no  two  adjacent  elements  overlap,  nor  are  there  any  gaps.  Also, 
that  any  one  of  four  types  of  EVD  elements,  (i.e..  Regular,  Leading  Edge,  Hinge, 
and  Infinity),  can  be  selected. 

Given  close  uniform  spacing,  the  program  would  always  yield  a good, 
smooth  chordwise  loading.  But  the  linear  nature  of  loading  given  by  the 
Regular  EVD  elements  allows  much  larger  spacing  to  be  used  in  regions  where 
the  loading  is  expected  to  nearly  linear,  that  is,  where  the  slope  of  the 
loading  is  nearly  constant.  In  addition,  in  order  to  more  accurately  approxi- 
mate the  chordwise  loading,  special  EVD  elements  have  been  developed.  These 
EVD  elements  provide  an  accurate  representation  of  chordwise  loading  in  these 
special  regions  which  are  normally  difficult  to  approximate  with  the  Regular 
EVD.  These  regions  are  near  the  leading  edge,  near  the  abrupt  turning  angle 
produced  at  a flap  or  jet  "hinge"  point,  and  far  downstream  on  the  jet  sheet. 
Use  of  the  special  EVD  elements  makes  chordwise  spacing  less  critical  and 
reduces  the  required  number  of  elements  on  each  section.  Most  of  the  chord- 
wise  EVD  elmenents  are  of  the  regular  type,  representing  a linear  loading 
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variation.  Each  section  will  contain  the  proper  arrangement  of  EVD  elements, 
automatically  assigned  by  the  program.  The  four  EVD  element  types  are  discussed 
in  detail  in  Volume  I of  this  report,  and  will  not  be  further  discussed  here. 

Chordwise  spacing  may  also  be  affected  by  section  camber.  In  regions 
where  there  is  a large  variation  in  camber,  close  spacing  should  be  used  in 
order  to  accurately  account  for  the  rapidly  changing  boundary  conditions 
resulting  from  the  wing  camber.  Volume  I of  this  report  contains  further 
discussion  and  examples  of  spacing  arrangements. 

Since  the  sectional  geometric  characteristics  of  a non-planar  wing  are 
approximated  by  the  "mean  camber"  line,  it  is  important  that  the  user  be  able 
to  translate  a typical  wing  section,  into  a suitable  camber  distribution.  See 
Figure  2 for  the  camber  line  representation  of  various  complex  wing  sections. 

All  of  the  wing  sections  originate  with  some  sfnooth  camber  distribution  of 
the  "clean"  configuration.  However,  when  various  types  of  control  surfaces 
are  deflected,  radical  changes  may  occur  in  the  shape  of  the  camber  line.  In 
each  case,  the  best  camber  representation  of  the  particular  deflected  configu- 
ration must  be  used.  Generally,  the  original  basic  camber  should  be  considered 
separately  from  the  changes  produced  by  deflections.  Ideally,  the  basic  camber 
distribution  would  be  fixed,  and  the  variations  in  camber  due  to  deflections 
would  be  superimposed  on  it.  The  program  cannot,  however,  automatically  take. 

# .f  * 

account  of  such  effects  as  ch^rd  extension  associated  with  flap  rotation  and 
slot  openings.  The  user-  should,  therefore,  include  some  allowance  for  those 
effects  in  the  basic  camber  distribution,  and  then  superimpose  s uple  deflec- 
tions which  do  not  contribute  to  changes  in  the  basic  camber.  Even  though 

the  program  may  be  used  to  approximate  such  complicated  configurations  as  ir> 

• ' 

Figure  2,  the  user  is  reminded  that  these  approximations  may  be  very  crude, 

*! 

and  rigorous  treatment  of  such  details  is  beyond  the  scope-  of  the  linearized 
approach.  Thus,  the  program  results  cannot  be  expected  to  accurately.predict 
the  aerodynamic  characteristics  of  such  configurations. 

2.4  Use  of  Fundamental  Cases 


An  important  feature  of  the  EVD  jet-wing  computer  program  is  the  concept 
of  fundamental  cases.  The  use  of  fundamental  cases  is  intended  to  increase  the 
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FIGURE  2.  REPRESENTATION  OF  WING  SECTIONS  BY  EQUIVALENT  CAMBER 
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amount  of  information  produced  by  the  program  and  to  make  the  use  of  that 
information  easier  for  the  user.  The  theory  of  fundamental  cases  is  discussed 
in  detail  in  Volume  I of  this  report,  and  will  not  be  repeated  here.  The  pro- 
gram application  and  treatment  of  fundamental  cases  will,  however,  be  explained. 

Fundamental  cases  are  used  to  treat  variations  in  different  types  of 
deflections  which  result  in  a change  of  wing  camber  and  are  not  applicable  to 
any  changes  in  planform  geometry.  The  types  of  deflections  which  may  be 
treated  as  fundamental  cases  are  as  follows: 

1 . Wing  twist 

2.  Jet  deflection 

3.  Camber 

4.  Leading  and  trailing  edge  flaps 

5.  Leading  edge  vertical  displacement  (special) 

The  fundamental  cases  produce  loadings  which  are  linearly  related  to 
the  magnitude  of  the  deflection  selected.  The  loading  of  each  case  may, 
therefore,  be  directly  added  to  the  loading  of  any  of  the  other  fundamental 
cases  to  produce  the  composite  loading  of  a particular  configuration.  The 
composite  cases  will  be  discussed  further  below. 

As  treated  by  the  program,  each  EVD  element  is  a part  of  every  funda- 
mental case.  Associated  with  each  element  is  a characteristic  set  of  incidence 
angles  relative  to  the  freestream  direction.  Each  of  these  incidence  angles 
is  the  resultant  angle  due  to  the  deflection  pattern  of  one  fundamental  case. 
When  the  deflection  of  a fundamental  case  varies,  the  associated  incidence 
angle  of  every  EVD  element  and  the  associated  loading  on  every  element  vary 
accordingly. 

An  example  of  the  use  of  fundamental  cases  for  a complex  wing  is  shown 
in  Figure  3.  In  Figure  3a,  the  planform  is  broken  up  into  regions  which  will 
be  deflected  as  control  surfaces.  The  program  allows  the  use  of  up  to  ten 
fundamental  cases,  the  first  of  which  is  automatically  set  up  as  a flat  plate 
wing  at  1°  angle  of  attack.  The  remaining  nine  fundamental  cases  may  be 
used  in  any  manner,  and  eight  of  them  have  arbitrarily  been  used  in  the  figure 
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b.  SEPARATE  FUNDAMENTAL  CASES  OF  A TYPICAL  WING  SECTION  (SECTION  A-A) 
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c.  FORMATION  OF  A TYPICAL  COMPOSITE  CASE  (SECTION  A A) 
FIGURE  3.  TREATMENT  OF  FUNDAMENTAL  AND  COMPOSITE  CASES 
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to  define  flaps,  jets,  leading  edge  flap,  camber,  etc. 

Figure  3b  shows  the  sectional  arrangement  of  the  same  fundamental  cases 
for  section  A-A.  It  will  be  noticed  that  in  each  case  the  incidence  of  all 
elements  outside  the  region  affected  is  zero.  Thus  when  the  deflection  is 
varied,  these  elements  will  have  no  change  in  incidence.  However,  since  the 
deflection  of  even  one  EVD  element  produces  a loading  distribution  on  every 
element  of  the  jet-wing  system,  variation  of  the  deflection  of  a fundamental 
case  will  produce  a corresponding  linear  variation  in  the  loading  of  every 
EVD  element. 

In  most  problems,  the  user  will  find  it  best  to  consider  those  components 
which  are  to  be  deflected  as  separate  fundamental  cases.  However,  it  is  some- 
times convenient  to  combine  several  types  of  deflection  into  one  fundamental 
case.  For  example,  camber  and  twist  will  usually  be  held  constant  while  other 
deflections  vary,  so  it  is  natural  to  combine  them.  However,  once  they  are 
combined,  they  can  no  longer  be  separately  distinguished  by  the  program. 

2.5  Use  of  Composite  Cases 

The  purpose  of  the  composite  case  capability  in  the  program  is  to  pro- 
vide the  user  with  data  for  any  configurations  that  may  be  required,  at  the 
minimum  computing  cost.  Figure  3c  shows  how  linear  fundamental  cases  are 
superimposed  in  any  manner  to  form  a composite  case.  The  result  is  the  com- 
posite loading  of  the  required  configuration.  The  composite  case  capability 
is  very  effectively  used  in  parametric  studies,  where  the  effects  of  systematic 
variation  of  parameters  are  required.  As  mentioned  in  the  discussion  of  fun- 
damental cases,  it  is  often  convenient  to  combine  several  unchanging  character- 
istics into  one  fundamental  case  which  may  then  be  included  in  all  composite 
cases  as  a "basic  wing"  starting  point.  The  user  supplies  the  scale  factors 
by  which  the  program  is  to  multiply  the  incidence  and  loading  of  each  respec- 
tive fundamental  case  before  adding  them  to  get  the  composite  incidence  and 
loading.  The  user  is  again  cautioned,  when  using  composite  cases,  about  ex- 
ceeding the  range  of  deflections  where  linear  behavior  is  likely.  He  should 
also  be  careful,  when  approximating  complex  sectional  geometric  characteristics 
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such  as  those  of  a slotted  flap  system,  to  limit  the  range  of  deflections  so 
as  not  to  significantly  change  the  planform  by  chord  extension. 

2.6  Final  Comments 


The  accuracy  of  the  results  obtained  with  this  computer  program  is 
generally  increased  as  the  number  of  EVD  elements  increases,  assuming  that 
the  elements  are  wisely  located  by  the  user.  When  pursuing  accuracy,  however, 
the  relative  significance  of  factors  which  are  unaccounted  for  must  be  kept 
in  mind.  In  addition,  the  computing  time  is  very  strongly  affected  by  the 
number  of  elements  (somewhere  between  a quadratic  and  cubic  function  of  the 
number  of  elements).  Thus,  a compromise  must  be  reached  between  the  desire 
for  the  best  program  results  on  the  one  hand  and,  on  the  other,  the  practical 
limitations  of  computing  resource:-  and  various  unpredictable  influences. 


It  is  again  strongly  recommended  that  the  user  read  this  computing 
manual  thoroughly  before  attempting  the  first  problem,  and  also  that  he  read 
Volume  I for  valuable  background  material  and  insight  into  the  theoretical 
foundations  on  which  the  program  is  based. 
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3.0  PROGRAM  DESCRIPTION 


3.1  General  Description 

The  McDonnell  Douglas  EVD  Jet-Wing  Lifting  Surface  Program  is  composed 
of  several  basic  component  programs,  each  of  which  is  independent  of  the  others 
in  logical  program  flow.  These  components  rely,  however,  on  each  other  for 
sharing  of  information,  and  all  of  them  are  required  for  satisfactory  analysis 
of  a problem.  The  component  programs  are  tied  together  by  two  master  routines, 
only  one  of  which  is  used  in  any  given  run.  The  choice  of  which  master  routine 
to  execute  is  made  by  the  main  routine  according  to  the  user's  requirements. 
Figure  4 shows  the  general  functional  arrangement  of  the  various  program  com- 
ponents. The  first  master  routine,  called  APPLY1 , controls  execution  of 
ordinary  runs,  without  stability  derivative  calculations.  If  stability  deri- 
vatives are  requested  by  the  user,  the  main  routine  calls  APPLY2,  which 
controls  execution  for  stability  derivative  runs.  Each  of  the  master  routines 
calls  subroutines  in  the  four  component  programs  as  they  are  needed. 

In  Component  1,  the  problem  is  formed  according  to  the  user's  input. 

All  input  is  read,  consistency  and  correctness  of  input  are  checked  as  far  as 
possible,  and  all  geometric  parameters  are  defined. 

In  Component  2,  the  mathematical  problem  is  solved.  The  system  of 
simultaneous  linear  equations  is  formed  using  the  downwash  influence  coeffi- 
cients of  each  EVD  element.  The  matrix  system  is  then  solved  directly,  and 
the  solution  is  achieved  in  the  form  of  the  vorticity  strength  of  each  EVD 
element  of  the  jet  wing  system. 

In  Component  3,  the  vorticity  solution  is  translated  into  useful  aero- 
dynamic coefficients  and/or  stability  derivatives.  These  include  chordwise 
loading,  spanwise  distribution  of  lift,  drag  and  pitching  moment,  and  total 
lift,  drag,  and  static  and  dynamic  moments. 

Component  4 contains  supplemental  subroutines  which  are  used  during 
calculation  of  dynamic  stability  derivatives.  This  is  not  a self-contained 
component  program  as  are  the  others,  but  is  simply  an  unstructured  collection 
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of  subroutines  whose  functions  are  related. 


3.2  Subroutine  Descriptions 

The  general  introduction  to  the  program  system  has  been  completed.  Next 
we  will  review  the  significant  details  and  functions  of  all  of  the  subroutines 
of  each  component  program.  Since  the  subroutines  are  generally  executed  in 
the  order  which  follows,  this  review  also  will  provide  a general  description 
of  the  logical  flow  during  execution.  Figure  5 gives  a very  general  flow 
chart  of  the  program  logic.  Reference  to  this  diagram  may  help  the  reader  to 
understand  the  overall  program  operation  and  the  function  of  the  various  sub- 
programs described  below. 

Main  Routine  (MAIN) 

The  main  routine  is  very  short,  and  serves  three  functions.  It  first 
reads  the  title  card  and  the  cards  containing  planform  parameters  and  control 
flags.  Next  it  checks  whether  stability  derivative's  have  been  requested  and 
accordingly  calls  either  APPLY!  or  APPLY2.  When  control  is  returned  to 
MAIN,  it  checks  whether  the  run  has  been  completed  normally,  whether  a new 
jet  strength  case  is  expected,  or  whether  a fatal  error  has  occurred. 

Master  Routine  for  Ordinary  Runs  (APPLY1) 

This  routine  is  very  short,  and  consists  of  simply  calling  the  control 
routines  for  component  programs  1,2,  and  3. 

Master  Routine  for  Stability  Derivative  Runs  (APPLY2) 

As  in  APPLY1,  the  control  routines  for  component  Programs  1,  2,  and  3 
are  called.  In  addition.  Component  4 is  called  to  produce  supplementary  infor- 
mation for  stability  derivatives  due  to  pitching.  The  loading  solutions  for 
the  first  run  are  then  stored  on  the  mass  storage  file  for  later  use,  and  new 
equivalent  camber  information  is  generated  for  each  element  due  to  yawing  and 
rolling  effects.  Component  2 is  again  called  for  solution  of  the  matrix  system. 
Finally  Component  3 is  called  for  integration  of  the  loading  to  get  dynamic 
stability  derivatives. 


16 


EVD 

ELEMENT 

GEOMETRY 


BOUNDARY 

CONDITIONS 


WING  PLANFORM 
AND  JET-SHEET 
GEOMETRY 


FUNDAMENTAL 

CASE 

GEOMETRY 


EVO 

ELEMENT 

GEOMETRY 


BOUNOARY 

CONDITIONS 


PLANFORM 

GEOMETRY 


FUNDAMENTAL 

CASE 

GEOMETRY 


J L. 


JET 

STRENGTH 


COMPOSITE 

CASE 

INPUT 


INPUT  (READ): 

• PLANFORM  GEOMETRY 

• JET  SHEET  GEOMETRY 

• FUNDAMENTAL  CASE 
GEOMETRY 

• COMPOSITE  CASE 
REQUIREMENTS 

• JET  STRENGTH 

OUTPUT: 

• EVO  ELEMENT  GEOMETRY 

• EVO  ELEMENT  TYPE 

• FUNDAMENTAL  CASE 
BOUNDARY  CONDITIONS 


PROBLEM  SOLUTION 


OOWNWASH 

INFLUENCE 

COEFFICIENTS 


H 

REGULAR  EVD 

L.E.  EVD 

U 

JET  INFINITY 
EVO 

HINGE  EVO 


REGULAR 


SINGULAR 


SINGULAR 

HINGE 

INFLUENCE 

COEFFICIENTS 


■ - 1 

r i 

r 

ICT 

WING 

WING 

JET 

* 

«*C  1 

INPUT  (INTERNAL): 

• EVO  ELEMENT  GEOMETRY 

• EVO  ELEMENT  TYPE 

• FUNDAMENTAL  CASE 
BOUNOARY  CONDITIONS 


OUTPUT. 

• VORTICITY  SOLUTION 
FOR  EACH  FUNDAMENTAL 
CASE 


COEFFICIENT 

SOLUTION  OF 

COLUMN 

HIM  1 HI  a 

N 

A.j'T'i.-B,. 

■i. 

STAGED!  ) AERODYNAMIC  QUANTITIES 


FUNDAMENTAL 

CASES 


CHOROWISE 

LOADING 


SPANWISE 

COEFFICIENTS 

ce-  Cm’  cd;-  C» 


COMPOSITE 

CASES 


STAGED! 

INPUT  (INTERNAL! 

• VORTICITY  SOLUTION 
FOR  EACH  FUNDAMENTAL 
CASE 

OUTPUT: 

AERODYNAMIC  COEFFICIENTS 

• CHOROWISE.  4C, 

• SPANWISE:  Ct.Cm.Cd 

• T0TAL:  Cl  Cm,CDi,<:,Cn,CJ.CY 


TOTAL  WING 
COEFFICIENTS 
CL'  CmCOi'Cl’  Cn'  CJ'CY 


PRINTED 

OUTPUT 


FIGURE  5a.  PROGRAM  STRUCTURE  FOR  A BASIC  RUN 


« t 


BASIC  RUN 


Control  of  Problem  Formation  (STAGE1) 


The  control  subroutine  (SGMAIN)  for  dealing  with  sectional  geometry  input 
and  element  definition  is  called  after  the  validity  of  certain  control  flags 
has  been  checked.  Alternate  methods  of  input  and  element  definition  may  be 
the  user-supplied  optional  routines  for  any  specified  values  of  the  IGTYPE 
flag.  (See  Section  4.1  - Input  Data  Instructions).  Next  the  composite  case 
requirements  are  read  by  calling  INCOMP.  Subroutine  BLOWIN  is  called  to 
read  the  jet  strength  requirements,  and  finally  BOXJ  is  called  to  prepare 
additional  jet  parameters  for  internal  use  later  in  the  program. 

Sectional  Input  Method  Control  (SGMAIN) 

The  spanwise  and  chordwise  spacing  of  wing  elements  are  read  by  calling 
subroutine  INPTS.  Next  the  planform  is  defined  by  calling  either  XLETR1 
for  irregular  planforms,  or  XLETR2  for  simple  trapezoidal  planforms.  The 
planform  parameters  are  normalized  by  half  the  input  value  of  wing  span  in 
subroutine  N0RM1 . Inputs  defining  the  spanwise  and  chordwise  spacing  of  ele- 
ments on  the  jet  are  read  by  calling  subroutine  INPUTJ.  Subroutine  BOXS 
next  combines  all  of  the  element  and  planform  input  data  acquired  up  to  this 
point  in  order  to  define  the  network  of  elements  required  by  the  program.  At 
this  point,  the  jet-wing  planform  arrangement  of  EVD  elements  is  completely 
defined,  but  it  remains  to  specify  the  pattern  of  deflections  required  by  the 
user.  This  is  accomplished  by  a sequence  of  calls  to  subroutines  INCASE, 
BEECEE,  and  0UT1  for  each  required  fundamental  case.  The  inputs  describing 
the  elements  to  be  deflected  and  the  amount  of  deflection  are  read  in  INCASE. 
BEECEE  interprets  this  data  and  derives  certain  parameters  to  be  used  later  in 
the  program.  0UT1  prints  the  resulting  parameters  for  each  element. 

If  shortcut  methods  can  be  devised  for  easier  input  of  all  or  part  of 
the  geometrical  data  required  for  a class  of  wings  of  particular  interest, 
some  of  the  routines  above  may  be  replaced,  modified  or  made  optional,  as  the 
user  requires.  However,  the  result  must  be  the  same  data  as  produced  by  the 
present  system  in  the  correct  form  for  use  by  remaining  components  of  the 
program. 
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Wing  Element  Input  (INPTS) 


This  subroutine  reads  the  spanwise  spacing  arrangement,  and  then  the 
number  and  chordwise  spacing  of  elements  on  the  wing. 

Irregular  Wing  Planform  Input  (XLETR1) 

This  subroutine  reads  the  leading  and  trailing  edge  input  coordinates 
to  define  the  arbitrary  planform,  and  linearly  interpolates  where  necessary 
to  define  the  leading  and  trailing  edges  at  any  sections  for  which  no  input 
is  given. 

Trapezoidal  Wing  Planform  Input  (XLETR2) 

This  subroutine  reads  the  general  planform  parameters  for  simple  sym- 
metrical trapezoidal  wings,  and  calculates  the  leading  and  trailing  edge 
coordinates  at  each  required  spanwise  section. 

Wing  Parameter  Normalization  (NORM!) 

This  very  short  subroutine  simply  normalizes  all  the  dimensional  plan- 
form  parameters  by  the  wing  half-span.  From  this  point  on,  the  program  deals 
entirely  with  the  scaled  wing,  where  the  wing  span  is  two  units. 

Jet  Element  Input  (INPUTJ) 

In  a manner  similar  to  that  used  in  INPTS,  the  number  and  chordwise 
spacing  of  elements  on  the  jet  are  read.  Note  that  the  spanwise  spacing  of 
jet  elements  is  the  same  as  for  the  wing,  but  that  some  or  all  of  these  jet 
sections  may  have  no  chordwise  elements,  thus  effectively  eliminating  the  jet 
at  those  sections. 

EVP  Element  Definition  (BOXS) 

This  subroutine  uses  the  input  data  from  the  above  subroutines  to  com- 
pute the  length,  width,  and  location  of  each  EVD  element  on  the  wing  and  jet. 
In  addition,  the  appropriate  EVD  type  is  chosen  for  each  element,  and  several 
other  parameters  are  defined  for  later  internal  use  by  the  program. 
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Leading  and  Trailing  Edge  Sweep  (TANS) 


This  subroutine  computes  the  sweep  angle  of  the  leading  and  trailing 
edges  at  each  section.  This  data  is  needed  later  for  calculation  of  side 
force. 

Fundamental  Case  Input  (INCASE) 

This  subroutine  reads,  for  each  fundamental  case,  the  flags  indicating 
the  types  of  deflections  to  be  used,  then  the  magnitude  of  each  particular 
type  of  deflection.  The  types  of  fundamental  case  deflections  read  are  twist 
angles,  leading  edge  vertical  displacement  (for  use  later  in  pitching  moment 
calculations),  jet  deflection  angles,  camber  angles,  and  hinge  location  and 
deflection  angles. 

Fundamental  Case  Boundary  Conditions  (BEECEE) 

From  the  fundamental  case  input,  this  subroutine  computes,  for  each 
fundamental  case,  the  resulting  incidence  angle  of  each  EVD  element  relative 
to  the  freestream.  This  is  done  by  summing  the  accumulative  deflections  due 
to  twist,  camber,  and  flap  deflection  from  the  leading  to  the  trailing  edges. 
The  effect  of  jet  deflection  is  also  included  on  the  first  jet  element  of  each 
section.  Thus  the  effects  of  all  deflections  are  superimposed  and  the  indi- 
vidual effect  of,  for  example,  deflection  of  a particular  flap  is  no  longer 
distinguishable  by  the  program.  Therefore,  if  the  user  wishes  to  be  able  to 
see  the  effects  of  that  flap  deflection  separately  from  those  of  the  other 
deflections,  he  must  make  use  of  a separate  fundamental  case,  where  only  the 
flap  deflection  is  present. 

Fundamental  Case  Output  (OUT!) 

This  subroutine  prints  out,  for  each  fundamental  case,  the  location, 
size,  resultant  incidence  angle,  hinge  deflection  angle,  if  any,  and  EVD  type 
of  each  EVD  element.  Parameters  of  spanwise  significance  are  also  printed, 
including  wing  chord,  leading  and  trailing  edge  coordinates,  leading  edge 
height,  twist  angle,  etc. 
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Composite  Case  Input  (INCOMP) 


This  subroutine  reads  the  user's  requirements  for  composite  cases,  in 
which  the  linear  fundamental  cases  may  be  superimposed  in  any  combination  to 
form  a deflection  pattern  of  particular  interest. 

Jet  Strength  Input  (BLOWIN) 

This  subroutine  simply  reads  the  sectional  value  of  jet  momentum  coeffi- 
cient, c^,  for  each  section  which  has  been  specified  to  have  jet  elements. 

Additional  Jet  Parameter  Definition  (BOXJ) 

This  subroutine  defines  an  additional  jet  strength  parameter  for  later 
internal  use,  and  checks  on  the  consistency  of  the  input  c^  data.  The  jet 
strength  from  the  previous  cy  case,  if  any,  is  also  saved  for  later  internal 
use. 

Control  of  Problem  Solution  (STAGE2) 

This  routine  simply  calls  the  two  overlay  segments  which  set  up  the 
system  of  simultaneous  linear  equations  and  solve  the  matrix  system. 

Development  of  Simultaneous  Equations  (STG2D) 

The  large  square  left-side  matrix  of  downwash  influence  coefficients  is 
computed  by  calling  DWNWSH.  Next  this  matrix  is  augmented  in  SHUFL1.  If 
previous  cp  cases  have  been  run,  these  steps  are  skipped  and  a simple  reaug- 
mentation is  performed  in  SHUFL2.  For  each  fundamental  case,  the  right  side 
column  of  constant  boundary  conditions  is  formed  in  C0LUM1.  Since  the  down- 
wash  influences  of  the  logarithmic  singular  hinge  EVD  elements  depend  only  on 
the  hinge  turning  angles  and  thus  are  known,  these  influences  are  computed  in 
HINGE  and  superimposed  directly  on  the  right  side  matrix  in  C0LUM2.  Sub- 
routine PREP  then  prepares  the  complete  linear  equation  system  for  solution 
by  assembling  the  left  and  right  side  matrix  data  together  on  tape  in  the  pro- 
per form  for  use  by  the  solution  subroutine. 
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Downwash  Influence  Coefficients  (DWNWSH) 

This  subroutine  calculates  the  downwash  on  every  element  due  to  the 
influence  of  all  the  other  elements.  Each  element  is  selected  in  turn,  and 
the  influence  of  all  elements  on  it  are  computed  one  at  a time  by  calling  the 
appropriate  downwash  influence  function:  EVD1 , EVD2,  or  EVD3  (hinge  EVD 
not  yet  considered).  If  the  wing  is  symmetric  or  anti -symmetric,  the  appro- 
priate influence  coefficients  are  superimposed  (added  or  subtracted)  to  reduce 
the  size  of  the  linear  system  to  half  the  number  of  equations.  The  array  of 
influence  coefficients  on  the  particular  element  under  consideration  forms  one 
row  of  the  coefficient  matrix,  and  these  rows  are  stored,  one  by  one  as  they 
are  computed,  on  the  mass  storage  device.  Unit  1.  Rows  corresponding  to  ele- 
ments on  the  wing  are  stored  first,  then  rows  corresponding  to  elements  on  the 
jet. 

Regular  Downwash  (EVD1) 

This  function  computes  the  downwash  at  any  point  in  the  wing  plane  due 

to  a triangular  distribution  of  vorticity  on  an  element  located  anywhere  on 

the  wing  (see  Figure  1). 

Leading  Edge  Downwash  (EVD2) 

This  function  computes  the  downwash  at  any  point  in  the  wing  plane  due 

to  a square  root  singular  distribution  of  vorticity  on  an  element  located  at 

the  leading  edge  of  the  wing  (see  Figure  1). 

Jet  Infinity  Downwash  (EVD3) 

This  function  computes  the  downwash  at  any  point  in  the  wing  plane  due 
to  a quadratic  decaying  distribution  of  vorticity  on  a trailing  jet  element. 
This  element  begins  several  chords  behind  the  wing  trailing  edge  and  extends 
downstream  to  infinity,  where  the  vorticity  decays  to  zero  (see  Figure  1). 
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Hinge  EVP  (EVD4) 


This  function  computes  the  downwash  at  any  point  in  the  wing  plane  due 
to  a logarithmic  distribution  of  vorticity  on  an  element  located  on  the  hinge 
line  of  a deflected  flap,  jet,  or  leading  edge  flap  (see  Figure  1). 

Initial  Matrix  Augmentation  (SHUFL1) 

This  subroutine  reads  the  matrix  rows  corresponding  to  elements  on  the 
jet,  subtracts  adjacent  rows  from  each  other,  and  further  modifies  certain 
elements  near  the  main  diagonal  due  to  the  nonlinear  influence  of  jet  strength. 
The  augmented  rows  are  then  replaced  on  mass  storage  Unit  1 immediately  behind 
the  original  downwash  matrix  rows. 

Additional  Matrix  Augmentation  (SHUFL2) 

This  subroutine  reads  the  augmented  matrix  rows  from  Unit  1,  remodifies 
certain  elements  near  the  main  diagonal,  then  restores  the  rows  in  their  origi- 
nal locations. 

Column  Matrix  Formation  (C0LUM1) 

This  subroutine  defines  the  right  side  column  matrix  of  constant  bound- 
ary conditions.  For  elements  on  the  wing,  the  values  are  the  total  incidence 
angles  relative  to  the  freestream,  in  radians.  For  the  first  element  of  each 
jet  section,  the  value  is  the  total  trailing  edge  incidence  angle  plus  the 
jet  turning  angle,  in  radians,  and  for  all  other  jet  elements,  the  values  are 
zero.  Each  column  corresponds  to  one  fundamental  case. 

Hinge  Downwash  Influence  (HINGE) 

For  each  element  on  the  wing  and  jet,  the  downwash  influence  of  each 
hinge  singularity  distribution  is  computed  by  calling  EVD4,  added  or  sub- 
tracted according  to  symmetry  requirements,  then  multiplied  by  the  appropriate 
hinge  turning  angle.  The  influence  due  to  all  such  hinges  are  summed  as  they 
are  calculated,  so  that  one  hinge  influence  factor  is  produced  for  each  element 
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on  the  wing  and  jet.  This  is  done  for  one  column  (fundamental  case)  at  a 
time. 

Column  Matrix  Hinge  Augmentation  (C0LUM2) 

This  subroutine  adds  to  the  right  side  column  matrix  row  corresponding 
to  each  element  on  the  wing  and  jet,  the  hinge  dov/nwash  influence  factor  com- 
puted in  HINGE.  For  rows  corresponding  to  elements  on  the  jet,  additional 
modifications  are  made  to  the  hinge  factors  before  adding  them  to  the  right 
side  matrix  in  order  to  account  for  the  influence  of  jet  strength  on  the  vorti- 
city  of  each  jet  hinge  element.  These  operations  are  done  for  one  column 
(fundamental  case)  at  a time. 

Matrix  Preparation  for  Solution  (PREP) 

The  form  in  which  the  matrix  solution  subroutine  expects  to  find  data 
is  row  by  row,  left  side  then  right  side  on  each  row,  stored  on  a sequential 
scratch  file.  Thus  the  augmented  left  side  rows  are  read  one  by  one  from  the 
mass  storage  Unit  1 into  a "transfer"  array,  the  elements  of  the  right  side 
matrix  are  defined  as  the  last  values  of  the  transfer  array,  and  the  entire 
array  is  written  as  one  record  onto  Scratch  Unit  2.  In  this  way  all  the  rows 
corresponding  to  elements  on  the  wing  and  then  the  jet  are  assembled  or 
"concatenated"  and  written  to  Unit  2 to  form  the  entire  matrix  system. 

Control  of  Matrix  Solution  (STG2S) 

This  control  routine  calls  the  matrix  solution  routine,  MATRIX,  re- 
trieves the  vorticity  solution,  and  calls  the  back  substitution  check  subrou- 
tine if  requested. 

Matrix  Solution  (MATRIX) 

This  subroutine  solves  the  matrix  system  directly  by  triangularization. 

It  is  able  to  solve  very  large  matrix  systems  using  a relatively  small  amount 
of  core  storage  because  it  deals  with  only  a portion  of  the  system  at  one  time, 
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the  remainder  being  stored  on  scratch  files. 


Scratch  File  Writing  (SAVE) 

Because  of  the  large  amount  of  writing  to  scratch  files  required  in  the 
matrix  solution  this  subroutine  is  used  to  write  whole  arrays  as  records  with- 
out reference  to  subscripts,  thus  significantly  increasing  the  solution  speed. 

Scratch  File  Reading  (GETT) 

As  with  subroutine  SAVE,  reading  whole  arrays  without  reference  to 
subscripts  significantly  reduces  the  matrix  solution  time. 

Matrix  Back  Substitution  (BAKSUB) 

As  a check  of  the  matrix  solution  accuracy,  the  user  may  wish  to  call 
for  a back  substitution  check.  In  a row  by  row  manner,  the  augmented  left 
side  matrix  is  read  from  mass  storage  Unit  1 and  the  corresponding  elements 
are  multiplied  and  summed  to  form  the  right  side  matrix.  This  right  side 
matrix  is  printed  for  reference  and  should  agree  with  the  boundary  conditions 
printed  earlier. 

Control  of  Aerodynamic  Coefficients  (STAGE3) 

This  routine  controls  computation  of  all  aerodynamic  coefficients,  in- 
cluding stability  derivatives,  for  both  regular  runs  and  stability  derivative 
runs.  If  stability  derivatives  are  not  requested  in  the  current  run,  chord- 
wise  loading,  spanwise  loading,  and  total  coefficients  are  computed  by  calling 
STG3FC  and  STG3FS,  respectively,  for  each  fundamental  case.  Then  a summary 
table  of  total  coefficients  for  all  fundamental  cases  is  printed  by  calling 
STG3FT.  For  each  composite  case  requested,  the  composite  loading  routine, 

STG3C,  computes  and  prints  spanwise  loading  and  total  aerodynamic  coefficients. 
If  stability  derivatives  are  requested,  subroutines  FUNDER  and  C0MDER  are 
called,  which  compute  stability  derivative  coefficients  for  all  fundamental  and 
composite  cases,  respectively. 


Chordwise  Loading  (STG3FC) 

The  loading  on  each  element  is  computed  from  the  vorticity  solution  pro- 
duced in  the  problem  solution  component  program.  If  requested,  the  loading  is 
printed,  along  with  detailed  expansions  of  the  loading  on  all  singular  leading 
edge  and  hinge  EVD  elements. 

Spanwise  and  Total  Loading  Control  (STG3FS) 

Spanwise  variation  of  lift,  pitching  moment  and  induced  drag  are  computed 
by  calling  SLOAD.  Total  induced  drag  is  also  computed  by  a momentum 
analysis  by  calling  subroutines  SLOADG  and  TREFTZ.  Total  lift,  pitching 
moment,  and  rolling  moment  are  computed  in  TLOAD,  and  total  induced  drag 
(by  pressure  integration)  and  yawing  moment  are  computed  by  calling  TLOADX. 

If  requested,  all  the  above  coefficients  are  then  printed.  This  sequence  is 
repeated  for  each  fundamental  case. 

Summary  Table  (STG3FT) 

If  requested,  subroutine  STG3FT  prints  a table  of  all  total  aerodynamic 
coefficients  for  all  fundamental  cases.  Since  no  coefficients  have  been  com- 
puted for  unused  fundamental  cases,  they  are  initialized  to  zero  before  printing. 

Composite  Case  Loading  Control  (STG3C) 

The  composite  chordwise  loading  is  first  computed  by  superimposing  the 
loading  from  the  required  fundamental  cases,  each  multiplied  by  its  respective 
input  scale  factor.  The  loading  is  also  printed,  if  requested.  For  composite 
cases,  both  the  loading  at  zero  angle  of  attack  and  the  linear  variation  with 
angle  of  attack  must  be  computed.  For  the  nonlinear  coefficients,  the  vari- 
ation with  angle  of  attack  becomes  quadratic,  and  a new  "cross-product"  term, 
linear  in  angle  of  attack,  is  also  required.  The  computation  of  all  spanwise 
and  total  coefficients  corresponding  to  the  zero,  linear,  quadratic,  and  cross- 
product  angle  of  attack  terms  are  computed  by  calling  the  various  loading 
utility  subroutines.  All  spanwise  and  total  aerodynamic  coefficients  are  then 
printed,  including  a summary  table  of  all  total  coefficients.  Finally,  sub- 
routine TABLE  is  called  to  give  a variation  of  total  lift,  drag,  and  moments 
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with  angle  of  attack.  The  entire  sequence  of  computations  above  are  repeated 
for  each  composite  case  requested. 

Initialized  Data  (BLKDTA) 

The  block  data  subroutine  defines  an  array  of  label  names  for  use  in 
labeling  the  summary  tables  of  total  aerodynamic  coefficients. 

Singular  EVP  Loading  Expansions  (EXPLE,  EXPHl,  EXPH21 

The  loading  at  five  points  on  the  EVD  element  are  computed  for  leading 
edge  singularities  and  the  leading  and  trailing  parts  of  hinge  singularities, 
respecti vely. 

Spanwise  Loading  (SLOAD) 

This  subroutine  integrates  the  chordwise  loading  at  each  section  to 
arrive  at  the  sectional  values  of  lift,  induced  drag  and  pitching  moment. 
Because  of  the  general  nature  of  the  computations,  the  loading  for  either 
fundamental  or  composite  cases  can  be  computed,  depending  on  the  element 
loading  given. 

Spanwise  Cross-Product  Loading  (SLOADX) 

The  sectional  cross-product  values  of  induced  drag  are  computed  by  chord- 
wise  integration  of  the  loading  of  each  element.  This  subroutine  is  utilized 
only  in  computation  of  composite  cases. 

Spanwise  Vorticity  (SLOADG) 

This  subroutine  integrates  the  chordwise  loading  at  each  section  to  get 
the  total  vorticity  of  the  jet  wing  system  for  use  in  calculation  of  induced 
drag  by  a momentum  analysis.  Integration  is  from  the  leading  edge  to  the 
trailing  edge  for  unblown  sections,  and  from  the  leading  edge  to  infinity 
for  sections  with  a jet. 
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Total  Loading  ( TLOAD ) 


This  subroutine  performs  spanwise  integration  of  the  sectional  values 
of  lift,  jet  strength,  and  pitching  moment  to  produce  total  lift,  induced 
drag,  jet  strength,  and  pitching  and  rolling  moments.  Only  total  coeffi- 
cients for  fundamental  cases  are  computed. 

Total  Loading  (TLOADO) 

This  subroutine  computes  total  lift  and  pitching  and  rolling  moments 
for  all  composite  cases  at  zero  angle  of  attack.  The  coefficients  are  com- 
puted by  summing  the  fundamental  case  coefficients,  each  multiplied  by  its 
respective  input  scale  factor. 

Total  Loading  (TLOADX) 

This  subroutine  computes  total  induced  drag  (both  momentum 
wing  pressure  integral  methods)  and  yawing  moment  coefficients  at  zero  angle 
of  attack  by  spanwise  integration  of  the  appropriate  sectional  data.  Both 
fundamental  and  composite  case  coefficients  may  be  computed. 

Trefftz  Plane  Downwash  (TREFTZ) 

This  subroutine  computes  the  induced  downwash  at  the  Trefftz  Plane 
(infinity)  due  to  all  the  loading  of  the  complete  jet-wing  system,  for  either 

fundamental  or  composite  cases.  These  data  are  used  in  the  momentum  induced 
drag  method. 

Composite  Case  Summary  (TABLE) 

This  subroutine  computes  and  prints  a table  of  the  variation  of  lift, 
induced  drag  (momentum  method),  and  pitching,  yawing,  and  rolling  moments  * 
with  angle  of  attack.  The  printout  is  given  for  each  composite  case. 

Fundamental  Case  Stability  Derivatives  (FUNDER) 

This  subroutine  controls  calculation  of  stability  derivatives  for  all 
fundamental  cases.  Subroutine  STG3FC  is  first  called  to  generate  the  chord- 
wise  loading  for  the  stability  fundamental  cases.  Then  both  chordwise  and 
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spanwise  integration  are  performed  by  calling  subroutines  SUMIT1  and  SUMIT2, 
to  yield  stability  derivatives  due  to  rolling  and  yawing  rates.  For  the  last 
fundamental  case  (due  to  rolling  rate)  certain  stability  derivatives  are  de- 
rived directly  from  chordwise  and  spanwise  loading  in  the  normal  manner  by 
calling  STG3FS. 

Composite  Case  Stability  Derivatives  (COMDER) 

This  subroutine  controls  computation  of  stability  derivatives  for  all 
composite  cases.  First  the  composite  case  loading  is  assembled,  requiring 
reading  of  the  fundamental  case  solutions  for  the  first  run  from  mass  storage 
Unit  1.  Then  the  subroutines  SUMIT1  and  SUMIT2  are  called  for  chordwise 
and  spanwise  integration  of  the  loading  to  obtain  the  required  stability  de- 
rivatives. Finally  a table  is  printed  showing  the  components  of  the  stability 
derivatives  and  the  variation  of  the  derivatives  with  angle  of  attack. 

Stability  Derivative  Integration  (SUMIT1) 

This  subroutine  computes  the  derivatives  of  yawing  moment  due  to  rolling, 
and  rolling  moment  due  to  yawing  by  chordwise  and  spanwise  integration  of  the 
loading  on  each  element. 

Stability  Derivative  Integration  (SUMIT2) 

This  subroutine  computes  the  derivative  of  yawing  moment  due  to  yawing 
by  chordwise  and  spanwise  integration  of  the  loading  on  each  element. 

Stability  Derivative  Table  (STABLE) 

This  subroutine  prints,  for  each  composite  case,  a summary  of  all  stabi- 
lity derivatives  and  the  terms  used  to  compute  each.  An  angle  of  attack  table 
is  then  computed  and  printed,  containing  all  the  stability  derivatives  which 
depend  on  angle  of  attack. 

Control  of  Utility  Routines  for  Stability  Runs  (STAGE4) 

During  the  first  run  of  a stability  derivative  run  sequence,  an  extra 
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fundamental  case  is  created  with  induced  camber  intended  to  simulate  a rate  of 
pitching.  This  is  accomplished  by  calling  subroutine  BCPICH,  and  all  the 
resulting  geometric  data  is  printed  by  calling  0UT2.  During  the  second  run 
of  a stability  derivative  run  sequence,  the  fundamental  case  solutions  from 
the  first  run  are  saved  on  mass  storage  Unit  1 in  subroutine  SAVECP.  The 
fundamental  case  geometry  is  then  redefined,  with  Induced  camber  simulating 
rates  of  yawing  and  rolling,  by  subroutines  BCYAW  and  BCROLL,  respectively. 
Again,  0UT2  gives  a record  of  the  new  fundamental  case  geometry. 

Fundamental  Case  Output  (0UT2) 

This  subroutine  is  the  same  as  subroutine  0UT2,  and  prints  the  geo- 
metric data  defining  the  stability  derivative  fundamental  cases. 

Save  First  Run  Solution  (SAVECP) 

For  stability  derivative  runs,  the  loading  solution  from  the  first  run 
must  be  saved  for  use  in  the  second  run.  This  is  accomplished  by  writing  the 
data  on  mass  storage  Unit  1,  directly  behind  the  matrix  information  previously 
stored. 

Pitching  Rate  Induced  Camber  (BCPICH) 

This  subroutine  defines  the  induced  camber  angles  on  all  wing  and  jet 
EVD  elements,  which  result  from  simulation  of  the  wing  pitching  about  the 
input  center  of  gravity  location. 

Rolling  Rate  Induced  Camber  (BCROLL) 

This  subroutine  defines  the  induced  camber  angles  on  all  wing  and  jet 
EVD  elements,  which  result  from  simulation  of  the  wing  rolling  about  the 
x-axis. 

Yawing  Rate  Induced  Camber  (BCYAW) 

This  subroutine  defines  the  induced  camber  angles  on  all  wing  and  jet 
EVD  elements,  which  result  from  simulation  of  the  wing  yawing  about  the  z-axis. 
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3.3  EVD  Functions-Numerical  Restrictions 


Because  of  the  characteristic  of  all  digital  computers  of  defining 
numbers  to  only  a limited  number  of  significant  digits,  it  sometimes  happens 
that  equations  cannot  be  evaluated  with  the  required  accuracy.  This  is 
particularly  a problem  where  the  difference  of  two  numbers  of  nearly  the 
same  magnitude  is  taken.  As  the  two  numbers  approach  each  other,  the  dif- 
ference loses  accuracy,  and  if  they  are  identical  (to  the  number  of  digits 


the  computer  can  hold),  the  difference  loses  all  accuracy. 


For  the  four  EVD  downwash  influence  functions  (see  Appendix  I of 
Reference  1),  the  above  type  of  problem  has  been  encountered  in  their  evaluation 
when  the  control  point  is  in  the  far- field  or,  for  certain  unique  locations, 
in  the  near-field  region  of  any  EVD  element.  These  can  generally  occur  along 
the  leading  and  trailing  edges  of  any  EVD  element,  and  at  the  apex  of  a 
Regular  EVD. 

The  following  restrictions  and  approximations  have  been  adopted,  without 
loss  of  generality,  in  order  to  prevent  random  inaccuracies  from  affecting  the 
results.  In  the  case  of  far-field  cutoff,  the  true  limiting  values  of  each 
function  are  used,  and  the  only  approximations  are  that  these  limiting  values 
are  used  somewhat  "closer"  than  real  infinity,  where  they  correctly  apply. 

For  near-field  cutoff,  the  limiting  values  are  used  in  small  regions  near  the 
points  where  the  values  correctly  apply.  The  subscript  i refers  to  the 
control  point  at  which  the  downwash  is  being  computed,  due  to  the  vorticity 
of  element  j. 
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Note:  Values  of  Hinge  EVD  regular  part  are  equivalent  to  Regular  EVD. 


4.0  INPUT  DATA  INSTRUCTIONS 


The  program  input  consists  of  various  types  of  cards  containing  all  the 
information  required  to  define  a particular  problem.  The  number  of  cards  of 
each  type  may  vary  from  problem  to  problem,  but  the  sequence  is  the  same  in  all 
runs.  As  a result  of  options  chosen  by  the  user,  some  types  of  cards  may  not 
be  required,  but  the  sequence  of  the  remaining  cards  must  not  be  changed.  Load 
sheets  for  all  cards  are  shown  in  Appendix  D.  The  input  deck  is  shown  diagram- 
matically  in  Figure  6. 

For  symmetric  and  anti -symmetric  wings,  all  input  and  output  data  are 
for  the  right  half  of  the  wing  only.  Figure  7 identifies  some  of  the  planfortn 
coordinates  required  as  input.  These  will  be  discussed  in  detail  below. 
Sectional  input  normally  begins  at  the  right  wing  tip  section,  working  inboard 
section  by  section.  However,  most  data  required  for  description  of  the  jet 
is  only  needed  at  sections  where  the  jet  exists,  and  unblown  sections  should 
be  skipped  as  indicated  in  the  instructions  below. 

Figure  8 shows  the  reference  coordinate  system  used  for  both  input  and 
output  data.  This  is  a wind  axis  system  with  the  freestream  direction  always 
aligned  with  the  x-axis.  It  must  be  noted  that  if  a wing  is  to  be  analyzed  in 
the  yawed  position,  it  must  be  input  as  a non-symmetric  planform,  rotated  about 
the  wing  apex  (z-axis)  with  respect  to  the  x-axis.  All  wing  sections  will 
still  be  aligned  with  the  x-axis,  including  the  wing  tip  sections.  The  printed 
output  will  also  be  referenced  to  the  wind  axis  system,  and  the  user  must  make 
suitable  transformations  if  coefficients  are  required  in  a body  axis  system. 


In  the  input  description  which  follows,  all  input  values  are  real  float- 
ing point  numbers  (F  format),  unless  otherwise  specified.  All  specified  inte- 
gers should  be  right- justified  in  their  fields. 
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4. 1 Input  Description 

• Title  Card  - This  card  provides  any  desired  description  of  the  computer 
run.  The  title  will  be  printed  at  the  top  of  the  first  page  of  output. 
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Column 


Name 


Explanation 

1-80  TITLE  Run  title  (any  acceptable  machine  characters). 

• General  Planform  Parameter  Card  - This  card  contains  basic  planform  infor- 


mation. 

Column 

Name 

Explanation 

1-10 

AREA 

Wing  area,  in  units  of  (SPAN)  to  be  used  for  normalization 

of  aerodynamic  coefficients. 

11-20 

SPAN 

Wing  span,  in  any  desired  units. 

21-30 

CREF 

Wing  reference  chord,  to  be  used  for  normalizing  various 
aerodynamic  coefficients.  Same  units  as  SPAN.  If  a value  of 
0.0  is  input,  the  mean  aerodynamic  chord  will  be  used. 

31-40 

XMC 

Point  about  which  pitching  moments  will  be  taken,  measured 
from  the  wing  apex  (origin).  Same  units  as  SPAN. 

41-50 

XCG 

Wing  center  of  gravity  location,  measured  from  wing  apex 
(origin),  which  will  be  used  as  pitching  axis  for  compu- 
tation of  stability  derivatives  due  to  pitching.  This  para- 
meter must  be  input  only  if  IDERIV  ? 0.  Same  units  as  SPAN. 

General 

Control 

Card  - This  card  contains  control  flaqs  which  describe  the 

basic 

characteristics  of  the  run. 

Column 

Name 

Explanation 

1-2 

NR0WS 

Number  of  spanwise  sections  (rows)  into  which  the  wing  is 

divided.  For  symmetric  or  anti -symmetric  wings,  only  the 
number  of  sections  on  the  right  half  should  be  input 
(3 s NR0WS  ^ 40).  Integer. 

3-4  NCASES  Total  number  of  fundamental  cases.  It  must  be  noted  that 
the  angle  of  attack  case  is  always  set  up  automatically  as 
fundamental  case  number  one,  and  no  input  data  is  required 
for  it.  Therefore,  NCASES  must  be  one  more  than  the  number 
of  cases  for  which  input  data  will  be  given,  to  allow  for 
the  angle  of  attack  case.  (1  £ NCASES  £ 10) . Integer. 

38 


- — - — — M 


. 


Column  Name 


Explanation 


| 


5-6 


7-8 


J 


10 


12 


14 


i 

1 


ISYMM  Symmetry  indicator  flag  (symmetry  about  the  x-axis). 

Integer. 

= 0 Wing  and  jet  are  symmetric. 

> 0 Wing  or  jet  are  non-symmetric. 

< 0 Wing  and  jet  are  anti -symmetric. 

I PRINT  Printed  output  control  flag.  Integer. 

> 1 Print  geometry  details  and  total  aerodynamic 

coefficients. 

= 1 In  addition,  print  spanwise  loading. 

= 0 In  addition,  print  chordwise  loading. 

< 0 In  addition,  print  all  matrices,  back  substitution 

check,  and  other  details.  This  option  is  normally 
reserved  for  troubleshooting,  since  it  produces  a 
very  large  amount  of  output. 

JETFLG  Jet  indicator  flag.  Integer. 

= 0 There  is  a jet  sheet. 

= 1 There  is  no  jet  sheet,  and  no  jet  inputs  will  be  read. 

IGTYPE  Wing  planform  geometry  indicator  flag.  Integer. 

= 1 Wing  planform  is  completely  arbitrary,  and  sectional 

leading  and  trailing  edge  coordinates  will  be  read 
to  define  planform. 

= 2 Wing  planform  is  trapezoidal,  and  simplified  plan- 

form  input  will  be  read. 

IHINGE  Hinge  EVD  indicator  flag.  Integer. 

= 0 Regular  EVD  only  will  be  used  on  all  hinge  elements. 

> 0 Hinge  EVD  will  be  used  on  all  hinge  elements.  This 

option  is  not  permitted  for  use  in  computing  dynamic 
stability  derivatives  (IDERIV>0). 
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Col umn  Name 


16 


IDERIV 


Explanation 

Dynamic  stability  derivative  control  flag.  Integer. 

= 0 A basic  run  will  be  executed,  with  no  stability 
derivatives  computed. 

> 0 A basic  run  will  be  executed,  and  in  addition, 

a dynamic  stability  derivative  run  will  be  executed. 
This  option  requires  approximately  double  the  basic 
computing  time. 


• Section  Centerline  Location  Cards  - These  cards  contain  the  spanwise  locations 

of  the  centerline  of  each  wing  (and  jet)section.  Eight 
values  per  card,  maximum  of  five  cards  (40  sections)  allowed. 


Column  Name 

1-10,11-20,  Y 
etc. 


Explanation 

Spanwise  distance  from  wing  centerline  (x-axis)  to  section 
centerline,  normalized  by  SPAN/2.  All  values  must  be 
(-1 .0  < Y < 1.0).  NR0WS  values  must  be  input,  beginning 
at  the  right  wing  tip  and  working  to  the  following: 

(a)  Wing  centerline,  for  symmetric  or  anti -symmetric 
wings. 

(b)  Left  wing  tip,  for  non-symmetric  wings. 


• Wing  Section  Type  Card  - This  card  indicates  the  chordwise  arrangement  of 

EVD  elements  for  each  section  on  the  wing. 


Column  Name 


1-2,  3-4, 
etc. 


ICTYPE 


Explanation 

Type  number  of  each  wing  section.  The  arrangement  of 
chordwise  EVD  elements  in  a row  (both  number  of  elements 
and  the  x/c  location  of  each)  determines  the  wing  row  type. 
Any  two  sections  which  have  the  same  number  and  x/c  loca- 
tions of  all  EVD  elements,  have  the  same  ICTYPE  value. 
NR0WS  values  must  be  input.  There  may  not  be  more  than  ten 
different  values  (i.e.,  section  types).  The  highest  value 
input  is  called  NWTYPE,  and  all  values  less  than  NWTYPE 
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Column  Name 


Explanation 


must  be  used,  that  is,  no  "gaps"  are  allowed  in  the 
sequence  1 through  NWTYPE.  Integer. 


• Number  of  Chordwise  Wing  Elements  - This  card  contains  the  number  of  chordwise 

EVD  elements  for  each  wing  section  type. 


Col umn 

Name 

1-2, 3-4, 

NI 

etc. 

Up  to 

CC10 

Explanation 

Number  of  chordwise  EVD  elements  of  each  wing  section  type. 
There  must  be  NWTYPE  values  (1  ^ NWTYPE s 10)  in  the 

sequence  for  section  type  1,  2,  3 (i.e.,  ascending 

order  of  row  types).  At  least  two  but  not  more  than  twenty 
chordwise  EVD  elements  must  be  input  on  each  wing  section. 
(2snI^20).  Integer. 


• Wing  Chordwise  Element  Coordinates  - These  cards  contain  the  x/c  coordinates 

of  each  EVD  element  of  each  section  type.  NWTYPE  sets 
of  cards  required,  each  with  NI  values  of  x/c.  Maximum 
of  ten  sets,  three  cards  per  set. 


Column  Name 


Explanation 


1-10,11-20,  XB  The  chordwise  coordinate  of  each  EVD  vortex  point,  measured 
etc.  from  the  sectional  leading  edge  and  normalized  by  the 

sectional  chord.  The  vortex  point  is  defined  as  the  leading 
edge  for  all  Leading  edge  EVD's  and  the  "peak"  point  for 
all  Regular  and  Hinge  EVD's.  The  first  value  of  each  set 
must  be  0.0  (leading  edge)  and  the  last  value  must  be  less 
than  1.0.  A maximum  of  20  values  (i.e.,  20  EVD  elements)  is 
permitted  for  each  set,  but  the  total  number  of  elements 
on  the  wing  and  jet  combined  must  not  exceed  600. 
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• Leading  and  Trailing  Edge  Coordinates  - These  cards  contain  the  coordinates 

of  the  leading  and  trailing  edges  of  the  given  sections. 
Only  stations  on  each  side  of  breaks  in  the  leading  or 
trailing  edges  need  be  input.  The  tip  and  root  sections 
must  be  input.  The  leading  and  trailing  edges  are  inter- 
polated for  sections  not  given  by  putting  straight  lines 
between  those  sections  v/hich  are  given.  This  method  of 
planform  definition  is  used  only  if  IGTYPE  = 1.  Number 
of  cards  required  is  2 2 and  sNRfJWS. 


Col umn 


Explanation 


1-10  Y Spanwise  distance  from  a section  centerline  to  the  center 

line  of  the  wing  (x-axis),  normalized  by  the  half  span. 
Each  value  must  be  the  same  as  one  of  those  already  input 
on  the  section  centerline  location  cards.  Values  on  the 
right  wing  half  are  positive,  and  values  on  the  left  wing 
half  (nonsymmetric  only)  are  negative. 


11-20  XLEAD  Chordwise  distance  from  section  leading  edge  (at  section 

centerline)  to  the  wing  apex  (y-axis).  Same  units  as  SPAN. 


21-30  XTRAIL  Chordwise  distance  from  section  trailing  edge  (at  section 

centerline)  to  the  wing  apex  (y-axis).  Same  units  as  SPAN. 

1 9 A 9 must  appear  in  CC1  of  the  next  card  after  all  the 

above  coordinate  cards  to  signal  that  all  desired 
sections  have  been  input.  This  card  is  required  only 
if  IGTYPE  = 1. 


• Trapezoidal  Wing  Parameters  - This  card  contains  planform  information  for 

simple  trapezoidal  wings.  It  replaces  the  above  coordi- 
nate cards  including  the  9 card  when  IGTYPE  = 2. 


Card  Variable 

Column  Name  Explanation 


1-10 

ARATI0 

Wing  aspect  ratio,  (SPAN'VAREA). 

11-20 

SWEEP 

Sweep  angle  of  wing  quarter-chord  line,  in  degrees. 

/ \ 

21-30 

TR 

Wing  taper  ratio,  CH0RDwing/cH0RDaxis  of  \ 

\ tip  / symmetry/ 

• Jet  Row  Type  Card  - 

This  card  indicates  the  chordwise  arrangement  of  EVD 
elements  for  each  section  on  the  jet  sheet.  Required 
only  if  JETFLG  = 0 . 

Column 

Name 

Explanation 

1-2, 3-4, 
etc. 

IJTYPE 

The  type  number  of  each  jet  section.  This  variable  is 
similar  to  the  ICTYPE  variable,  except  that  for  sections 
with  no  jet  a value  of  0 must  be  input.  The  number  of 
different  jet  section  types  is  called  NJTYPE.  The  number 
of  nonzero  values  input  is  NROWSJ,  the  number  of  sections 
having  a jet.  Integer. 

• Number 

of  Chordwise 

Jet  Elements  - This  card  contains  the  number  of  chord- 

wise  EVD  elements  for  each  jet  section  type.  Required 
only  if  JETFLG  = 0 . 

Column 

Name 

Explanation 

1-2,  3-4 
etc. 

, NI 

Number  of  chordwise  EVD  elements  of  each  jet  section 
type.  This  variable  is  similar  to  NI  for  wing  sections 

above,  except  that  there  must  be  NJTYPE  sections.  At 
least  2 but  not  more  than  10  chordwise  EVD  elements  must 
be  input  on  each  jet  section  (2«NIS10).  Integer. 

• Jet  Chordwise  Element  Coordinates  - These  cards  contain  the  x/c  coordinates 

of  each  element  of  each  jet  section  type.  NJTYPE  sets 
of  cards  required,  each  with  NI  values  of  x/c.  Maximum 
of  10  sets,  2 cards  per  set.  Required  only  if  JETFLG  = 0. 
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Column  Name 


Explanation 


1-10,11-20,  XB  The  chordwise  coordinate  of  each  EVD  vortex  point, 

etc*  measured  from  the  sectional  leading  edge  (at  centerline) 

and  normalized  by  the  sectional  chord.  The  first  value 
of  each  set  must  be  1.0  (trailing  edge)  and  a maximum 
of  10  values  is  permitted  for  each  set  (2<NI<10). 

• Fundamental  Case  Control  Card  - This  card  identifies  the  types  of  linear 

geometric  variations  to  be  included  in  each  fundamental 
case.  The  number  of  fundamental  cases  input  must  be 
one  less  than  NCASES,  to  allow  for  the  angle  of  attack 
case.  A new  fundamental  case  control  card  is  required 
for  each  of  the  (NCASES  - 1)  input  cases.  In  each  of 
the  flags  below,  a zero  value  indicates  omission  of  the 
respective  type  of  input  for  that  fundamental  case.  A 
non-zero  value  indicates  the  variation  will  be  included 
and  input  must  be  given  to  define  it. 


Col umn  Name 


2 

4 

6 

8 

10 


INTWST 

INHITE 

INDELJ 

INCAMB 

INBETA 


Explanation 

Spanwise  twist  distribution  flag.  Integer. 

Leading  edge  vertical  displacement  flag.  Integer. 
Jet  deflection  flag.  Integer. 

Camber  flag.  Integer. 

Wing  hinge  deflection  flag.  Integer. 


• Twist  Distribution  Card  - These  cards  contain  the  spanwise  distribution  of 

wing  twist.  NR0WS  values  required,  eight  per  card. 
Required  only  if  INTWST  f 0. 


Col umn  Name 


1-10,11-20  TWIST 
etc. 


Explanation 

Sectional  wing  twist,  in  degrees,  at  section  centerline, 
with  respect  to  the  wing  reference  plane.  Positive  is 
in  the  same  sense  as  a positive  angle  of  attack.  See 
Figure  3b  (2). 
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• Leading  Edge  Height  Card  - These  cards  contain  the  displacement  of  the  wing 

leading  edge  from  the  wing  reference  plane.  NR0WS  values 
required,  eight  per  card.  Required  only  if  INHITE  t 0. 


Col umn 


1-10,11-20,  HO 
etc. 


Explanation 

Displacement  coordinate  of  sectional  leading  edge  from  the 
wing  reference  plane,  normalized  by  the  sectional  chord. 
These  values  are  used  only  for  computation  of  the  moment 
arm  of  the  jet  reaction  thrust  contribution  to  pitching 
moments.  Leading  edge  displacement  may  be  the  result  of 
dihedral,  twist,  nonlinear  movement  of  a leading  edge 
device,  etc.  Translation  resulting  from  ordinary  linear 
leading  and  trailing  flap  deflections  and  angle  of  attack 
are  accounted  for  automatically  by  the  program. 


• Jet  Deflection  Card  - These  cards  contain  the  spanwise  variation  of  jet 

deflection,  relative  to  the  trailing  edge.  NR0WSJ 
values  are  required,  eight  per  card.  Required  only  if 
INDELJ  f 0. 


Column 


Explanation 


1-10,11-20,  DJ 
etc. 


Jet  turning  angle,  in  degrees,  relative  to  the  mean  line 
of  the  trailing  edge.  Positive  deflection  is  downv/ard. 
See  Figure  3b  (D • 


• Camber  Type  Flag  Card  - This  card  indicates  the  chordwise  distribution  of 

camber  for  each  section  on  the  wing.  Required  only  if 
INCAMB  ? 0 . 


Column  Name 

1-2, 3-4  ICT 
etc. 


Explanation 


Camber  type  number  of  each  wing  section.  The  arrangement 


..  f 


of  camber.  Including  all  positions  (x/c)  and  angles  must 
be  the  same  in  order  for  two  sections  to  have  the  same 
value  of  ICT.  NR0WS  values  must  be  input,  but  a maxi- 
mum of  10  different  types  is  allowed.  The  highest  value 
input  is  called  NCT,  and  there  may  be  no  "gaps"  in  the 
numbering  sequence  between  0 and  NCT.  A zero  value 
indicates  no  camber.  Integer. 

I Camber  Angle  Cards  - These  cards  contain  the  camber  angles  for  each  camber 

section  type.  NCT  sets  of  cards  required,  three  cards 
per  set  maximum.  Required  only  if  INCAMB  f 0. 

Explanation 

The  camber  angle,  in  degrees,  at  the  downwash  control 
point  of  each  EVD  element  ( i . e. , the  point  arbitrarily 
chosen  as  halfway  between  any  two  adjacent  XB  points, 
including  the  trailing  edge).  Positive  values  are  in 
the  same  sense  as  positive  angles  of  attack  See  Figure 
3b  <D. 

• Trailing  Edge  Camber  Angle  Card  - These  cards  contain  the  trailing  edge  angle 

due  to  camber  only.  NR0WSJ  values  are  required.  Required 
only  if  both  INCAMB  i 0 and  JETFLG  = 0. 

Explanation 

Trailing  edge  angle,  in  degrees,  due  to  camber  only.  These 
values  are  required  only  for  sections  which  have  a trail- 
ing jet,  and  are  used  only  in  computing  the  total  jet  de- 
flection angle  with  respect  to  the  freestream. 

• Hinge  Section  Type  Card  - This  card  identifies  the  arrangement  of  leading 

and  trailing  flap  hinges  on  each  section.  NR0WS  values 
are  required.  Required  only  if  INBETA  f 0. 


Column 

Name 

1-10,11-20 

ACTE 

etc. 

Column  Name 

1-10,11-20,  AC 
etc. 
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Column 


1 -2,3-4  IHT 
etc. 


Explanation 

Hinge  section  type  flag.  If  there  are  no  hinges,  a value 
of  0 should  be  input.  Sections  with  the  same  hinge  sec- 
tion type  must  be  alike  in  number  of  hinges,  their  loca- 
tion (x/c),  their  type  (leading  or  trailing  flap)  and  in 
all  deflections.  There  may  be  as  many  as  NR0WS  different 
section  types.  The  number  of  different  types  is  called 
NHT,  and  there  may  be  no  gaps  in  the  sequence  0 through 
NHT.  Integer. 


• Hinge  Parameter  Cards  - These  cards  contain  the  location,  type  (leading  or 

trailing  edge  flap)  and  turning  angle  at  each  hinge  on  a 
given  type  of  hinge  section.  NHT  cards  required.  Each 
section  may  have  a maximum  of  four  hinges.  Required  only 
if  INBETA  1 0. 


Column  Name 

1-10,21-30  XBH 
etc. 


11,31, 

etc. 


12-20,32-40,  BETA 
etc. 


Explanation 

Distance  from  the  sectional  leading  edge  to  the  hinge 
point,  normalized  by  the  sectional  chord.  This  value 
must  be  the  same  as  one  of  the  XB  coordinates  input  for 
the  section. 

Hinge  type  identification  flag.  Integer. 

= 0 Trailing  flap  hinge  (positive  deflection  in  the 
sense  of  positive  angle  of  attack). 
f 0 Leading  flap  hinge  (positive  deflection  in  the 
sense  of  negative  angle  of  attack). 

Hinge  deflection  angle,  in  degrees,  of  the  element 
behind  the  hinge  point,  relative  to  the  element  before 
the  hinge  point.  See  Figure  3b  ®,  ©,  ®. 


Composite  Case  Cards  - These  cards  specify  the  desired  superposition  of  the 

linear  fundamental  cases.  A maximum  of  24  composite  cases 
may  be  requested.  No  composite  cases  may  also  be  specified 
(9  card  alone  requirea).  See  Figure  3c. 
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Column 


Explanation 


Name 


1-2,9-10,  N 

etc. 

3-8,11-16,  A 

etc. 

1 9 


Fundamental  case  number  which  is  to  be  included  in 
forming  a given  composite  case.  Integer. 

Multiplication  factor  to  be  applied  to  the  fundamental 
case  identified  in  the  previous  two  card  columns. 

A "9"  must  appear  in  CC1  of  the  card  following  the  com- 
posite case  cards  to  signal  the  end  of  composite  case 
input.  (Required  even  if  there  are  no  composite  case 
cards  input). 


• Jet  Strength  Cards  - These  cards  contain  the  jet  strength  for  all  sections 

which  have  a jet.  NROWSJ  values  are  required.  An 
unlimited  number  of  sets  of  cards  may  be  input,  maxi- 
mum of  five  cards  per  set.  Required  only  if  JETFLG  = 0 . 


Column 


Name 


Explanation 


1-10,11-20,  CMU 
etc. 


Sectional  jet  momentum  strength  for  each  jet  row.  CMU 
is  defined  as  CMU  = J/qc(y),  where  J is  the  sectional 
jet  momentum  per  unit  span  and  q and  c are  the  dynamic 
pressure  and  sectional  chord,  respectively.  There  may  be 
no  zero  values  input,  unless  all  values  are  zero  (indi- 
cating no  jet  at  all). 

A "9"  must  appear  in  CC1  of  the  card  following  all  CMU 
cards  to  signal  the  end  of  CMU  input.  Required  only  if 
JETFLG  = 0 . 


4.2  Input  Restrictions 

Most  of  the  limitations  imposed  on  the  range  of  input  values  are  designed 
to  minimize  computer  storage  requirements,  but  still  to  provide  the  capability 
for  adequate  treatment  of  complicated  configurations.  While  all  of  the  input 
restrictions  have  been  stated  or  implied  in  the  input  description  section  above, 
they  will  be  summarized  here  for  convenient  reference. 
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1 . 3 £ NR0WS  £ 40 

2.  1 £ NCASES  £ 10 

3.  IGTYPE  = 1 or  2 

4.  Number  of  chordwise  wing  elements  on  any  section  2 £ NI  £ 20 

5.  Number  of  chordwise  jet  elements  on  any  section  2 £ NI  £ 10 

6.  Maximum  of  ten  different  section  types  each  on  wing  and  jet 

7.  On  wing  0.0  £ XB  < 1.0 

8.  On  jet  1.0  £ XB 

9.  Only  NROWSJ  values  required  for  DJ,  ACTE  and  CMU  input 

10.  0.0  £ CMU  < 800.0 

11.  For  stability  derivative  runs,  only  one  set  of  CMU  data  may  be  input 


' 

ij 

i 


Due  to  certain  curve  fitting  restrictions  in  the  calculation  of  induced 
drag  by  the  momentum  method,  at  least  three  adjacent  sections  must  be  used 
in  defining  either  a jet  sheet  or  a region  on  the  wing  which  has  no  trailing 
jet  sheet.  Thus,  a jet  cannot  consist  of  only  one  or  two  sections.  Likewise, 
inboard  or  outboard  of  any  jet,  there  must  be  at  least  three  unblown  wing 
sections.  If  the  three-in-a-row  rule  is  not  followed,  the  program  will  stop 
automatically. 


I c 

! 

I t 
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5.0  OUTPUT  DATA  DESCRIPTION 


The  program  output  consists  entirely  of  data  written  to  the  standard 
print  Unit  6.  Four  print  options,  controlled  by  the  input  value  of  the  IPRINT 
flag  described  above,  are  available  to  the  user.  The  normal  printout  is  given 
for  the  default  value  of  zero,  and  includes  chordwise  and  spanwise  loading  and 
total  wing  aerodynamic  coefficients.  If  the  details  of  chordwise  or  spanwise 
loading  are  not  required,  a large  part  of  the  print  and  a slight  amount  of  com- 
puting time  may  be  saved  by  requesting  one  of  the  more  restricted  forms  of  output. 

If  more  information  is  required  for  program  checkout  or  troubleshooting, 
the  extended  print  option  may  be  used.  Under  this  option,  all  of  the  perti- 
nent data  of  the  problem  solution  are  printed,  including  downwash  coefficients, 
augmented  matrix  coefficients,  vorticity  solution  and  back  substitution  check 
for  recovery  of  the  right  side  matrix.  The  matrix  information  is  given  in  a 
row  by  row  sequence,  each  row  corresponding  to  one  EVD  element.  For  each  row 
the  influence  coefficients  of  all  EVD  elements  on  the  one  in  question  are 
listed  across  the  page.  For  example,  the  downwash  influence  due  to  Element  12 
on  Element  6 would  be  the  second  number  of  the  second  line  listed  under  matrix 
Row  6. 


The  detailed  print  option,  when  used  in  a stability  derivative  run,  will 
also  give  for  the  second  run  the  chordwise  loading,  spanwise  loading  of  the 
last  fundamental  case,  and  the  summary  table  of  total  coefficients  for  all 
fundamental  cases.  These  data  are  not  intended  for  direct  use,  but  mainly 
for  indirect  checking  of  the  consistency  of  certain  stability  derivative 
results. 

It  will  be  helpful  for  the  user  to  refer  to  the  sample  problem  output 
of  Appendix  A during  discussion  of  the  output  data  given  below.  All  output 
symbols  and  labels  are  listed  in  the  nomenclature  section  and  all  angles 
are  measured  In  degrees. 

Even  though  many  significant  digits  are  printed,  this  is  by  no  means 
meant  to  Imply  that  the  coefficients  are  accurate  to  such  a degree.  Six  or 
eight  significant  digits  are  printed  only  so  that  the  effects  due  to  small 
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deflections  or  planform  changes  will  be  visible  in  the  F format.  While 
changes  in  the  last  significant  digits  may  indicate  the  trends  of  effects  due 
to  small  changes,  they  must  not  be  construed  to  indicate  their  magnitude  ac- 
curately. 

The  composite  case  linear  aerodynamic  characteristics  are  lift,  pitching 
moment,  root  bending  moment,  and  rolling  moment.  They  are  composed  of  a value 
at  zero  angle  of  attack  and  a linear  variation  with  angle  of  attack,  and  may 
be  calculated  as  follows: 


C.  = C,  + C.  a 
L Q=0  La 


(CCL  = CCLO  + CCLA  * ALPHA) 


The  nonlinear  characteristics  include  all  of  the  drag  contributions,  yawing 
moment  and  side  force,  and  are  composed  of  a value  at  zero  angle  of  attack, 
a square  variation  with  angle  of  attack,  and  a "cross-product"  linear  vari- 
ation with  angle  of  attack  as  follows: 


CD  = CD. 


a=0 


+ CD  a + cr 

ua 


cr 


(CCD  = CCDO  + CCDX  * ALPHA  + CCDA2  * ALPHA  **  2) 


The  proper  methods  of  composition  of  the  dynamic  stability  derivatives  vary 
greatly,  and  are  shown  in  the  printout  for  each  of  the  types  of  derivatives. 


As  noted  in  the  input  data  instructions,  the  reference  system  of  both 
input  and  output  data  is  a wind  axis  system,  with  the  x-axis  aligned  with  the 
freestream.  If  it  is  required  to  analyze  a yawed  wing,  the  wing  must  be  input 
in  the  yawed  position  and  considered  as  non- symmetric.  If  aerodynamic  coeffi- 
cients are  required  in  a body  axis  system,  they  must  be  transformed  by  the 
user.  In  the  linearized  analysis,  the  only  contribution  to  side  force  should 
be  the  leading  edge  suction,  taken  normal  to  the  local  leading  edge. 


The  program  automatically  normalizes  all  the  geometry  data  by  dividing 
by  half  the  input  value  of  the  wing  span  (i.e.  SPAN).  Thus,  the  program  works 
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exclusively  with  a wing  which  has  a span  of  two  units,  and  all  output  is  for 
this  normalized  wing. 

5.1  Ming  Characteristics  and  Control  Flags 

The  program  title  header  and  the  input  case  title  are  printed  for  ref- 
erence. Next  the  major  wing  planfortn  characteristics  and  input  control  flags 
are  printed.  Since  the  program  will  automatically  redefine  some  of  the  para- 
meters and  flags  under  certain  run  conditions,  both  the  input  and  internally 
used  values  are  printed  for  reference. 

5.2  Fundamental  Case  Geometry 

The  page  header  is  printed  at  the  beginning  of  data  for  each  fundamental 
case.  For  each  section,  the  pertinent  dimensions  are  given  on  the  first  line. 
Next  the  detailed  data  for  the  wing  are  given,  with  the  general  information 
printed  on  the  first  line,  followed  by  the  chordwise  arrangement  of  all  ele- 
ment locations,  angles  and  EVD  types.  Each  number  printed  across  the  page 
corresponds  to  one  EVD  element,  beginning  at  the  sectional  leading  edge  and 
working  toward  the  trailing  edge.  The  data  given  include,  for  each  element, 
the  following: 

a.  Chordwise  location,  first  normalized  by  the  local  chord,  then 
in  dimensional  units. 

b.  Chordwise  element  length. 

c.  Total  incidence  angle  resulting  from  all  forms  of  fundamental 

case  deflections. 

d.  Hinge  turning  angle  of  each  element. 

e.  EVD  type. 

Next, the  jet  element  data  is  given,  if  there  is  a jet  at  that  section.  The 
first  line  contains  data  applying  to  the  whole  jet  section,  and  is  followed 
by  the  chordwise  arrangement  of  EVD  element  data.  The  element  incidence  angle 
is  not  printed  since  it  has  no  physical  meaning,  and  the  hinge  turning  angles 
must  be  zero  for  all  but  the  leading  jet  element.  The  length  of  the  jet  in- 
finity EVD  element  will  always  be  printed  as  a large  number  overflowing  the 
format,  since  it  Is  theoretically  infinite  and  Is  automatically  assigned  the 
value  1010. 
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The  values  of  TYPE 


identify  the  following  EVD  types: 


10 

20 

30 

41 

42 

43 


Regular  EVD 
Leading  Edge  EVD 
Infinity  EVD 

Leading  edge  flap  hinge  ] 

Trailing  edge  flap  hinge  > Hinge  EVD 
Trailing  edge  - jet- emersion  ] 


For  all  fundamental  cases  after  the  first  the  printout  of  element 
dimensional  locations  and  lengths  are  skipped,  since  they  are  the  same  for 
all  cases.  The  normalized  chordwise  locations  are  also  repetitious,  but  are 
printed  for  reference  in  each  case. 


5.3  Sectional  Jet  Strength 

This  output  is  given  for  reference  and  includes  the  c^  value  at  every 
section,  whether  or  not  it  is  a jet  section.  This  printout  is  skipped  if 
JETFLG  f 0. 


For  dynamic  stability  derivative  runs,  an  extra  fundamental  case  is 
created  automatically  by  the  program,  with  incidence  angles  defined  according 
to  the  quasi-steady  induced  camber  due  to  a one  degree/second  pitching  rate 
about  the  XCG  location. 


5.4  Fundamental  Case  Coefficients 

After  the  problem  solution  is  completed,  the  loading,  ACp,  on  each  EVD 
element  is  printed  for  all  fundamental  cases.  These  data  are  given  section 
by  section,  with  the  wing  element  loading  first  and  the  jet  element  loading 
second.  In  addition,  detailed  loading  at  small  intervals  are  given  for  all 
leading  edge  and  hinge  EVD  elements.  This  is  intended  primarily  as  an  aid 
to  the  user  in  plotting  the  chordwise  loading.  The  hinge  detailed  print  is 
given  only  if  IHINGE  i 0,  since  hinge  loadings  are  otherwise  represented  by 
regular  EVD  elements  rather  than  singular  hinge  EVD  elements. 


The  spanwise  loading  data  include,  for  each  fundamental  case,  the  lift, 


induced  drag,  pitching  moment,  and  center  of  lift  for  each  wing  section.  At 
the  bottom  of  each  column,  the  values  labeled  "TOTAL"  are  the  corresponding 
integrated  aerodynamic  coefficients  for  the  whole  wing.  The  sectional  pitch- 
ing moment  coefficients  are  taken  about  the  sectional  leading  edges,  while 
the  total  coefficients  are  taken  about  both  the  wing  apex  (origin)  and  about 
the  input  XMC  location.  Sectional  center  of  pressure  and  center  of  lift 
data  are  all  normalized  by  the  sectional  chord.  The  total  values  indicate  the 
chordwise  centers  of  pressure  and  lift,  each  normalized  by  either  reference 
chord  or  wing  semi  span.  The  total  value  of  induced  drag  computed  by  the 
momentum  method  is  given  under  the  column  headed  ALFINF. 

For  dynamic  stability  derivative  runs,  the  derivatives  due  to  pitching 
are  derived  from  the  last  fundamental  case  and  are  printed  below  the  spanwise 
loading  data. 

Finally,  a summary  table  of  all  total  aerodynamic  coefficients  for  all 
fundamental  cases  is  printed.  These  values  are  the  same  as  those  printed 
under  the  spanwise  loading  data,  but  also  include  rolling  and  yawing  moments. 
The  most  important  coefficients  are  labeled  with  asterisks. 

It  should  be  noted  that  in  order  to  compute  the  sectional  pitching 
moment  about  the  leading  edge  due  to  the  thrust  component  (in  the  freestream 
direction)  of  jet  reaction  acting  at  the  trailing  edge,  the  vertical  distance 
between  the  leading  and  trailing  edges  must  be  known.  For  wings  with  flap 
or  slat  fundamental  cases,  the  program  computes  this  distance  using  the  hinge 
points  and  respective  deflection  angles.  If,  however,  the  same  deflections 
are  described  by  a camber  fundamental  case,  the  program  has  no  way  of  knowing 
that  the  leading  or  trailing  edges  are  displaced,  and  therefore,  the  jet  thrust 
pitching  moment  contribution  will  always  be  zero.  The  user  is  thus  warned  to 
use  care  when  interpreting  the  pitching  moment  results  for  jet-wings  where 
flap  or  slat  deflections  are  treated  as  camber  alone. 

5.5  Composite  Case  Coefficients 

The  composite  case  output  is  similar  to  that  for  fundamental  cases,  ex- 
cept that  all  data  is  given  for  each  composite  case  before  going  on  to  the 


next.  The  first  line  of  data  contains  the  input  multiplication  factors  for 
the  fundamental  cases  comprising  this  particular  composite  case.  Then  the 
chordwise  loading  is  printed  for  each  section.  The  ACp  values  are  for  zero 
angle  of  attack,  labeled  (A  = 0),  and  for  the  variation  of  loading  with  angle 
of  attack,  labeled  (A  = 1).  The  fundamental  case  factors  apply  only  to  the 
loading  at  zero  angle  of  attack,  and  the  loading  variation  with  angle  of 
attack  is  taken  directly  from  fundamental  case  number  one. 

As  noted  above,  the  CP  value  is  infinite  at  the  Leading  Edge  and 
Hinge  points,  and  the  values  printed  are  the  Leading  Edge  "average  value" 
and  the  Regular  EVD  "peak"  value  underlying  the  singular  Hinge  EVD  loading, 
if  any.  Thus,  these  points  must  not  be  plotted  directly  by  the  user.  The 
Jet-Infinity  EVD  decays  slowly  from  its  "peak"  value  to  zero  at  infinity 
downstream.  The  printed  values  may  be  used  to  compute  any  intermediate 
loading  from  the  following  equations: 


Leading  Edge  Element: 

A c 


'XB 


,2  f/XB\*  _ /XB\1  + /XB\ 

3 pi  |\«i / \6i/J  pi+l\V 


Hinge  Element: 

H 


AC. 


4Bi 


"PXB 


Ac 


PXB 


log  c£.  , XB-XH 
1 og ( c [ XB-XH  | ) + ^-4^ ' 

Vi 

(for  XB  < XH) 

46i  log  (cS.)(XB-XH) 

1 og(c | XB-XH | ) ' 

(for  XB  < XH) 


AC, 


Ki  tcPH)(xB-XH) 

i i-l 


AC 


p,-  ^^(xb-xh) 


Jet-Infinity  Element: 


Ac  = Ac  d2  c(XB-l) 
PXB  pi  I V 'I 


n 

I 

where  c is  the  local  chord,  l is  the  element  length,  XB  is  the  point 
where  the  loading  is  required,  XII  is  the  hinge  point,  and  i refers  to 
the  Leading- Edge,  Hinge  or  Jet- Infinity  elements,  respectively.  XB,  XH, 
and  6 are  normalized  by  c. 

The  spanwise  loading  data  are  presented  in  generally  the  same  way  as 
for  each  fundamental  case,  except  that  there  is  more  information.  For  the 
linear  aerodynamic  coefficients  (lift,  pitching  moment,  and  center  of  lift), 
the  sectional  and  total  coefficients  are  given  for  both  zero  angle  of  attack 
(labels  ending  with  0)  and  their  variation  with  angle  of  attack  (labels 
ending  with  A).  The  first  two  lines  of  total  lift  center  data  are  the  chord- 
wise  locations  normalized  by  reference  chord,  and  in  the  third  and  fourth 
lines  they  are  normalized  by  the  wing  semi  span. 

For  the  nonlinear  induced  drag  coefficients,  three  terms  are  required: 
zero  angle  of  attack  (labels  ending  with  0),  variation  with  angle  of  attack 
squared  (labels  ending  with  A2),  and  the  "cross-product"  variation  with  angle 
of  attack  (labels  ending  with  X). 

A summary  table  of  all  total  aerodynamic  coefficients  is  printed,  with 
the  data  taken  from  the  total  data  of  the  spanwise  loading.  Coefficients  at 
zero  angle  of  attack  and  those  which  are  associated  with  a and  a2  are 
shown. 

Finally,  a summary  table  of  the  variations  of  total  lift,  pitching 
moment  (about  XMC),  rolling  moment,  induced  drag  (momentum  method)  and  yaw- 
ing moment  about  XMC  with  angle  of  attack  are  given.  In  the  sample  case,  the 
small  negative  drag  values  are  not  errors,  but  are  simply  an  indication  of 
inaccuracy  resulting  from  the  very  coarse  spacing  used  In  the  example. 

| 

* 

It  must  be  noted  that  in  the  above  description  of  composite  cases, 
angle  of  attack  equal  to  zero  refers  to  the  basic  configuration  of  the  wing 
as  constructed  from  the  scale  factors  input  for  all  fundamental  cases.  Since 
the  a- case  is  a legitimate  fundamental  case  and  a multiplication  factor  may 
be  input  for  Case  1,  it  is  possible  to  have  a composite  case  which  is  labeled 
"«a0"  by  the  program,  but  which  actually  has  an  arbitrary  constant  flat  plate 
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loading  superimposed  on  it.  The  program  does  not  distinguish  an  input  for 
fundamental  Case  1 from  the  other  cases  and,  therefore,  does  not  "know"  that 
the  wing  is  already  at  some  angle  of  attack. 

This  feature  may  be  used  with  care  to  obtain  the  composite  loading  at 
some  angle  of  attack  of  particular  interest,  say  a=4.3°,  but  the  user  should 
be  aware  that  he  has  now  forced  the  loading  labels  to  be  misleading.  A label 
of  a=0°  now  refers  to  the  loading  at  a=4.3°,  a=10°  now  means  a=14.3°, 
etc.  The  user  should  check  the  scale  factor  for  fundamental  Case  1 printed 
at  the  top  of  the  page  in  order  to  avoid  misinterpretation  of  the  composite 
case  output.  If  it  is  0.0,  then  all  the  labels  are  consistent  with  the 
printed  data. 

All  of  the  output  data  described  in  Section  5.5  are  repeated  for  each 
requested  composite  case  in  turn. 

5.6  Stability  Derivative  Output 

For  most  runs  the  program  will  repeat  the  above  output  for  all  requested 
composite  cases.  It  will  then  attempt  to  read  new  CMU  data  and,  if  no  CMU 
data  is  found,  it  will  attempt  to  begin  a completely  new  run.  If  no  new  run 
data  is  found  the  job  is  terminated.  For  dynamic  stability  derivative  runs, 
however,  the  program  will  automatically  begin  a second  run  in  which  the  de- 
rivatives due  to  rolling  and  yawing  are  computed.  All  previously  input 
fundamental  cases  are  altered  to  represent  the  effects  induced  by  a yawing 
rate  of  one  degree/second.  Likewise*  the  extra  fundamental  case  is  altered 
to  represent  the  effects  induced  by  a rolling  rate  of  one  degree/second. 

The  fundamental  case  output  for  the  second  run  consists  simply  of  the 
rolling  and  yawing  derivatives  due  to  yawing  (original  fundamental  cases)  and 
due  to  rolling  (extra  fundamental  case). 

The  output  for  the  composite  cases  again  includes  the  input  fundamental 
case  multiplication  factors.  As  with  the  basic  aerodynamic  coefficients, 
the  stability  derivatives  are  generally  composed  of  several  angle  of  attack 
terms,  but  some  of  them  are  also  dependent  on  the  rates  of  yawing  or  rolling 


themselves.  The  complete  formulas  for  calculation  are  given  In  the  printout. 

Finally,  a table  of  the  variation  of  those  derivatives  which  depend  on 
angle  of  attack  is  given.  Since  only  one  CMU  case  can  be  treated,  the  program 
will  attempt  to  begin  a completely  new  run  after  the  second  stability  deri- 
vative run.  The  termination  message  is  printed,  indicating  the  current  run 
is  completed,  and  if  no  further  input  data  is  given,  the  job  is  terminated. 


6.0  OPERATIONAL  CONSIDERATIONS 


6.1  Program  Versions  Available 

The  EVD  Jet-Wing  Computer  Program  has  been  written  entirely  in  Fortran 
IV,  and  could  run  on  any  large  computer  system.  The  program  was  developed  on 
the  IBM  360  and  370  systems  at  Douglas  Aircraft  Company  and  converted  for  use 
on  CDC  6000  series  computers.  It  is  expected  that  only  minor  changes  would 
be  required  for  running  on  any  other  large-scale  computing  system. 

6.2  Overlay  Structure 

Because  of  the  large  size  of  the  program,  the  overlay  feature  of 
FORTRAN  is  required  for  operation  on  most  machines.  Shown  below  in  Figure  9 
is  the  overlay  structure  of  the  Mark  II  version.  This  structure  requires 
approximately  141000  octal  words  of  storage  on  the  CDC  6500  and  approximately 
235000  bytes  of  storage  on  the  IBM  360  computer.  The  overlay  control  cards 
for  CDC  6000  series  computers  are  included  in  the  program  listing.  Appendix  C. 

6. 3 Peripheral  Storage  Devices 

The  program  requires  the  use  of  six  peripheral  storage  devices,  plus 
any  devices  required  by  the  local  system.  These  devices  may  be  tape  or  disk 
units,  and  are  described  as  follows: 


Unit  Type 


Use 


1 

2 

3 

4 

5 

6 


Binary  Mass  storage  of  downwash  influence  coefficients  and 

matrices. 

Binary  Matrix  input  to  solution  routine. 

Binary  Temporary  storage  during  matrix  solution. 

Binary  Temporary  storage  during  matrix  solution. 

BCD  Standard  input  unit. 

BCD  Standard  printed  output  unit. 
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6.4  Program  Running  Time 

The  computer  running  time  for  the  EVD  Jet-Wing  Program  depends  primarily 
on  the  number  of  EVD  elements  which  are  used  to  represent  the  jet-wing  and, 
of  course,  on  which  computer  is  being  used.  The  number  of  fundamental  and 
composite  cases  requested  has  a relatively  small  influence  on  computing  time. 

Figure  10  shows  the  time  required  for  program  execution  on  the  IBM 
370/165  computer.  CPU  (Central  Processor  Unit)  time  refers  to  real  time 
used  in  calculation.  10  (lnput/0utput)  time  refers  to  an  equivalent  time 
charged  for  processing  on  peripheral  storage  devices  and  may  be  knov/n  as  PP 
(Peripheral  Processor)  time.  This  includes  reading  and  writing  of  temporary 
data  sets  such  as  the  card  input  unit,  printed  output  unit,  scratch  tapes, 
punched  cards,  etc.  Total  time  is  the  sum  of  CPU  time  and  10  time. 

Other  machines  and  other  installations  may  require  a different  amount 
of  time  for  execution.  It  is  expected  that  CPU  time  will  vary  among  machines, 
but  will  be  representative  of  the  speed  of  any  particular  model.  On  the  other 
hand,  10  time  is  expected  to  vary  from  installation  to  installation,  depending 
on  the  facilities  available,  the  demand  placed  upon  them,  their  cost,  etc. 

Figure  10  applies  to  the  execution  of  one  value  of  jet  strength  (which 
may  be  the  no-jet  case).  For  each  additional  jet  strength  case  required  in 
a single  run,  the  total  time  will  be  increased  by  approximately  half  the  total 
time  of  the  first  jet  strength  case. 
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APPENDIX  A - SAMPLE  PROBLEM 


Both  the  input  card  Images  and  the  program  printout  for  a small  sample 
case  are  given  below.  It  must  be  noted  that  this  case  has  been  prepared  only 
as  a brief  example  of  program  input  and  output.  The  spacing  is  very  rough, 
both  spanwise  and  chordwise,  and  should  not  be  used  as  a model  for  preparation 
of  other  cases.  The  fundamental  and  composite  case  capabilities  of  the  pro- 
gram are  far  greater  than  those  demonstrated  by  the  sample  case.  Because  of 
the  crude  spacing  and  few  elements,  it  is  expected  that  the  results  shown 
will  not  accurately  predict  the  wing's  aerodynamic  characteristics. 
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SAMPLE  PROBLEM  OUTPUT 


• EVO  JET  - WING  COMPUTER  PROGRAM  • 


•••  ONR  SAMPLE  CASE  •*•  R PC  TAKGULAX  MING  CMU  ■ l WITH  STABILITY  DERIVATIVES 


AREA  ■ 
SPAN  » 
CREf  • 
XMC  ■ 
CMAC  • 
ARATIO  • 
XCG  • 


USED 

0. 888080 
2.000000 
0. 444444 
0. 11111! 
0.444444 
4.500000 
0.111111 


1 NPUT 
4.500000 
4.500000 
1 .000000 
0.250000 
0.999999 
*.500000 
0.250000 


NRCWS  - 4 

NCASES  ■ 3 

ISVMM  ■ 0 

1PR INT  • 0 

JETFLG  • 0 

JGTYPE  » 2 

IHINGE  • 0 


NUMBER  OP  WING  ELEMENTS  » 21 

NUMBER  OF  JET  ELEMENTS  • 16 

TOTAL  NUMBER  OP  ELEMENTS  - 37 


• ELEMENT  GEOMETRY  DATA  AND  FUNDAMENTAL  CASE  DATA  FOR  FUNDAMENTAL  CASE  1 • 


SECTION  1 Y * 3.97*000  DELTA  - 0. 025000  XLE AD 


XTRAIL  ■ 3.444444  CHORD  ■ 0.444444  TANLE  - 0.0 


WING  ELEMENTS 


0.0 

0.0 

0. IU0)30 
1.000000 
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TWIST  « 0.0 

0.100000  0.200000 
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0. I 00000 
1.000000 
0.0 


0.088889 
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0.0 

10 


HL  ■ 0.0 

0.500000 
0.222222 
0.400000 
1.000000 
0.0 
10 


THETA  S ■ 0.0 

0.900000 
0.4  00000 
0.100000 
1.000000 
0.0 


JET 

ELEMENTS 

NJ  « 4 

D « U. 888889  OJ 

a 0.0 

ACTE  - 

0.0  THETA 

■ 1.000000 

XH 

1 . Ou 0000 

1 .099999 

1.500000 

3.000000 

XI 

3 . 4444  44 

0.4BBP8C 

0. 066667 

1.333333 

TEL 

5.099999 

0.400 001 

1.500000* 

etTA 

5.0 

0.0 

0.0 

0.0 

TYPE 

10 

10 

10 

30 

... 

SECTION  2 •• 

* V a 0.8975U0  DELTA 

* 0.062500  XLE AO  « 

0.0 

XTRAIL  > 0.444444 

CHORD  - 0.444444 

TANLE  • 

0.0 

WING 

ELFMENTS 

NW  « 5 

TwIST  * 0. 

0 

HL  « 0.0 

THETA  S - O.J 

Xd 

0.0 

0. 100 000 

0.200000 

0.500000 

0.900000 

XI 
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0.0 

0.0 

0.0 
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T YPt 

20 

10 

10 

10 

10 

JET 

ELEMENTS 

NJ  * 4 

n * 0 ,°PRRP9  CJ 

* 0.0 

ACTE  - 

0.0  THETA 

« l. OOOOOO 

XB 

1 .050000 

1 .w999°«; 

1.  500000 

3 . 000000 

X I 

3.444444 

0.439^89 

0. 666667 

1.333333 

CEL 

3 .0Soc 90 

0.40000 l 

1. 500000* 

RET  A 

3.0 

0.0 

0.0 

0.0 

TYPE 

10 

10 

10 

30 

SECT ICN  3 *' 

* V s o. 

687500  DELTA 

* 0.1375 

00  XLE  AO  ■ 

0.0 

XTRAIL  ■ 0.444444 

CHORD  ■ 0.444444 

TANLE  • 

0.0 

Ml  NO 

ELEMENTS 

f;w  * 6 

TwIST  = 0. 

0 

H.  - 0.0 

THETA  S ■ 0.3 

XH 

3.0 

0.100000 

0.200000 

0.500000 

0.800000 

0.900000 

XI 

3.0 

0.044  '.4* 

0. 38*80° 
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0.355555 

0.400000 

DPI 

J.loOOJO 
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0.300000 

0.300000 
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l.OOOOOC 
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BETA 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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20 

10 

10 

10 

10 

10 

JET 
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0.0 

0.0 

0.0 

TYPE 

10 

10 

10 

30 

• ** 
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0.0 
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0.0 

WING 
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Nw  • 5 
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0.100000 
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0.500000 
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X 1 
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0.400000 
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0.100000 
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• ELEMENT  GE34ETRV  OAT  A ANO  FUNDAMENTAL,  CASE  OATA  FOR  FUNDAMENTAL  CASE  2 • 


•••  SECTION  I •••  T * 0.97*000  DEL  TA  ■ 0.02*000  ALE  AO  • 

MING  ELEMENTS  NM  • * TWIST  - 0.0  HL  - 0.0 


JET  ELEMENTS 


XB 

3.0 

0. 100000 

0.200000 

0.500000 

0. 900000 

EPS 

0.0 

0.0 

0.0 

0.0 

0.0 

BETA 

0.0 

0.0 

0.0 

0.0 

0.0 

TYPE 

20 

10 

10 

10 

10 

J ■ 4 0 ■ 0.888889  DJ  > 1.000000 

1.000000  1.099994  1. *00000  1.000000 

1.000000  0.0  0.0  0.0 

43  10  10  30 


XT  RAIL  • 0.4444*4  CHORD  ■ 0.4 

THETA  S • 0.0 


THETA  • 1.000000 


TANLE  ■ 0.0 


•••  SECTION  2 ••• 

MING  ELEMENTS  Nl 

Xft 
EPS 
BETA 
TYPE 

JET  ELEMENTS  N 

XB 
BETA 
TYPE 


Y • 0.887500  DELTA  ■ 0.062500  X LEAD  » 

X • 5 TWIST  • 0.0  HL  • 0.0 

0.0  0.100000  0.200000  0.500000 

3.0  0.0  0.0  0.0 

0.0  0.0  0.0  0.0 

20  10  10  10 

J » 4 0 » 0 ,888099  DJ  « 1.000000 

l.OOOJOO  1.099999  1.500000  3.000000 

1.0J03J0  0.0  0.0  0.0 

43  10  10  30 


XTRAIL  • J. 4444*4  CHORD  - 0.444444  TANLE  • 0.0 


THETA  S ■ 0.3 

0.400000 
0.0 
0.0 
10 


THETA  > 1.000000 


•••  SECTION  3 Y ■ 

Ml«  ELEMENTS  K»  • 6 

X6  0.0 
EPS  3.0 
BETA  3.0 
TYPE  20 


3.687503  CclTA 


0.137500  XLEAP 


XTRAIL  • 0.444444  CHORO  - 0.444444  TANLE  • 0.0 


TWIST  • 0.0  HL  ■ 0.0 

0.100000  0.200000  0.500000 

0.0  0.0  0.0 

0.0  0.0  0.0 

10  10  10 


THETA  S ■ 0.3 

o.eooooo  o.9ooooo 


JET  ELEMENTS 

xp 

BETA 

TYPE 

•••  SECTION  * «• 
MING  ELEMENTS 


J - 4 0 • 0 .088869  OJ  * 1.000000 

1.000300  1.099999  1.500000  3.300000 

1.000000  0.0  0.0  0.0 

43  13  10  30 

Y - 3.27*000  DELTA  - 0.275000  XLEAD  - 


THETA  ■ l. 000000 


XTRAIL  - 0.*44444  CHORD  • 0.444444  TANLE  ■ 0.0 

THETA  S • 0.3 


X3 

3.0 

0.130  GOO 

0.  200000 

0.500000 

0.900000 

EPS 

3.0 

0.0 

0.0 

0.0 

0.0 

BETA 

3.0 

0.0 

0.0 

0.0 

0.0 

TVPF 

20 

10 

10 

10 

10 

TS 

NJ  - 4 

0 • 0 .808089 

OJ 

« 1.000000 

ACTE  ■ 

XB 

L.OJOJOO 

1 .099999 

1.500000 

3.000000 

BETA 

1.000000 

0.0 

0.0 

0.0 

TVPF 

43 

10 

10 

30 

THETA  - 1.300000 


• ELEMENT  GEOMETRY  DATA  AND  FUNDAMENTAL  CASE  OATA  FOR  FUNDAMENTAL  CASE  3 • 


SECTION  1 •••  Y ■ 0.97*000  CEL TA  ■ 0.025000  XLEAD  • 

MING  ELEMENTS  nm  ■ 5 TwIST  - 0.0  HL  > 0.0 


1 

XB 

EPS 

3.0 

3.0 

0.100000 

0.0 

0.200000 

0.3 

3.500000 

0.0 

0.9000 00 
C.O 

* 

BETA 

3.0 

0.3 

0.0 

0.0 

0.0 

TYPE 

23 

10 

10 

10 

10 

XTRAIL  > 3.444444  CHORD  ■ 0.444444  TANLE  ■ 0.0 

THETA  S ■ 0.0 


JET  ELEMENTS  NJ  « * 

Xft  1.003300 
BiTA  3.0 
TYPE  *.3 


0 ■ 0.886EP9  DJ  » 0.0 

l.3999<c  1.500000  3.000000 

0.0  0.0  0.0 

10  10  30 


SECTION  2 Y - 0. 907*00  CELTA  » 0.062500  XLEAD  - 


nING  ELEMENTS  NM  ■ 5 

Xft  3.0 
EPS  3.0 
BETA  3.0 
TYPE  20 


TMIST  « 0.0  HL  - 0.0 

0. 100C00  0.200000  0.500000 

0.0  0.0  0.0 

U.J  0.0  3.0 

10  10  10 


■ 3.444444  CHORO  > 0.444444  TANLE  • 0.0 


THETA  S • 0.3 

0.900000 
0.0 
0.0 
10 


X*  l.OoOCOO  1.049999  1.500000  3.000000 

BETA  3.0  0.0  0.0  3.0 

TYPc  43  10  10  30 

•••  SECT ICN  3 • *•  V • 3.687*00  DELTA  » 0.137530  XLE AO  ■ 

M|»A  ELEMENTS  Nm  • 6 TwIST  * C.O  HL  > 0.0 


XB 

0.0 

0.  100000 

0.200000 

0.500000 

0.800000 

0.900000 

EPS 

3.0 

0.0 

0.0 

0.0 

0.0 

1 . 030000 

beta 

0.0 

0.0 

0.0 

0.0 

0.0 

1.000000 

TYPE 

20 

10 

13 

10 

10 

42 

JET 

ELEMENTS 

NJ  > 4 

o « o.see*89 

OJ 

• 0.0 

ACTE  ■ 

0.0 

>8 

1.000.* 

30  1.099999 

1.500000 

3.000000 

beta 

3.0 

0.0 

0.3 

0.0 

TYPE 

43 

10 

10 

10 

•••  SECTION  4 •• 

• V « 

J. 275003  DELTA 

» 0. 

,27*000  XLEAD  ■ 

0.0 

XTRAIL  • 

MING 

ELEMENTS 

NM  ■ 5 

TMIST  > 0. 

0 

HL  • 0.0 

THETA  S ■ 0, 

KE- 

3.0 

O.lOGCOO 

0. 200000 

0.5  00000 

0.900000 

EPS 

3.0 

0.3 

0.0 

0.0 

0.0 

BETA 

3.0 

3.3 

0.0 

•3.0 

0.0 

TYPE 

20 

10 

10 

10 

10 

JET 

ELEMENTS 

NJ  « 4 

C « 0.88*889 

DJ 

« 0.0 

ACTE  • 

0.0 

xe 

1.GJ0300  1.099*99 

1.500000 

1.000000 

BETA 

3.0 

0.3 

0.9 

0.0 

TVPF 

43 

10 

13 

10 

XTRAIL  ■ 3 . 444444  CHORD  ■ 0.444444  TANLE  • 0.0 

THETA  S ■ 0.3 


THETA  - 1.000000 


XTRAIL  « 3.444444  CH0*O  ■ 0. *44444  TANlE 


TMETA  • 0.0 


THIS  PAGE  IS  BEST  QUALIIT'  PRACTICABLE 
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• SECTIONAL  JET  PLOWING  COEFFICIENTS  • 


CHU 

1.000000 

1.000000 

1.000000 

1.000000 


* ELEMENT  GfOMETBY  data  AND  PUNOAMFNTAL  case  data  for  FUNDAMENTAL  case  a * 


•••  SECTION  1 


0.975000  TELTA  * 0.02*000  XLFAD 


WING  ElFMfMS 


S NW  * 5 TWIST  « 0.0 

xB  0.0  0.100000  0.200000 

EPS  -O.hOOJOO  -0.200000  0.199999 


BETA  0.0 
TYPE  20 

JET  ELEMENTS  NJ  » * 

Xh  1.000000 
BETA  0.0 
TYPE  43 


HL  « 0.0 

0.500000 
0.900000 
0.0 
10 


1.099*9°  1.500000  3.000000 


XTRAIL  > 3.44*4**  CHORD  ■ 0. ******  TANLE  ■ 0.0 


THETA  S - 0.3 

0.900000 
1.400000 
0.0 
10 


THETA  ■ 1.503000 


•••  SECTION  2 **♦  V « J.8975C0  DELTA 


0.062500  XL?  AD 


WING  ELEMENTS 


S NW  * 5 UIST  * 0.0  HL  * 0.0 

XF  3.0  0.100000  0.200000  0.500000 

EPS  -3.400000  -0.200000  0.199999  0.900000 


BIT/  3.0 
TYPI  23 

JET  ELEiMFMTS  NJ  « a 

XP  1.000000 
BET/  3.0 
TYPE.  43 


C « 0.RP8R89  DJ  * 0.0 

1.0O9499  1.503000  3.000000 

0.0  0.0  0.0 

10  10  30 


*••  SECTION  3 **'  Y » 0.69  7 500  DEL TA 


0.137500  XLFAD 


a I NO  ELEMENTS  NW  * 
XB  0.0 


TWIST  * 0.0 


.100000  0.200000  0.500000 


EPS  -3.400l.00  -0.200000  0.  199999 

BITA  3.0  0.0  0.0 

TYPE  20  10  10 


JET  ELEMENTS  NJ  ■ 4 

X B 1 .000,  00 
BETA  J.G 
TYPE  43 


D » U . 9B8BP9  DJ 

1.099094  1.  500000 

0.0  0.0 

10  10 


0.8  00000 

0.0 


0.0 

3.000000 

0.0 

30 


XTRAIL  ■ 3.444444  CHORO  ■ 0.444444  TANLE  ■ 0.0 


THETA  S - 0.0 

0.900000 
1.400000 
0.0 
10 


THETA  * 1.500000 


XTRAIL  ■ 3.444444  CHORD  • 0.444444  TANLE  » 0.0 


THETA  S ■ 0.0 

0.800000  0.900000 

1.199999  1.400000 

0.0  0.0 

10  42 


THFTA  ■ 1.503000 


WING  t LE  Mf  NTS 


»*»  Y * 0.27*000  PEL  TA  * 0.275000  XLFAO  ■ 


TwIST  * 0.0  HL  * 0.0 

0.  100  )00  0.200000  0.500000 


EPS  -0.-G3C00  -3.2JC000  0.  199999 

BETA  3.0  0.0  0.0 

TYPE  20  10  10 


JET  ELF  M(  NTS  NJ  ■ A 

X't  l.OuOJOO 
BETA  0.0 
TYPE  43 


0 « 0.838389  DJ  * 0.0 

1.099999  1.500000  3.000000 

0.0  0.0  0.0 

10  10  30 


XTRAIL  ■ 0.4*4444  CHORO  « 0.444444  TANLE  • 0.0 


THETA  S ■ 0.3 

0.900000 
1.400000 
0.0 
10 


THFTA  • 1.503000 
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* CtOaOMlSC  LOAOING  FOR  «U  FUNO«MfNTAl  CASES  • 


SECTION 

1 V • 

0.9  75000 

CHORD  > 

0.444444 

1 

xb 

C*Sl  1 

CASE  2 

CASE  3 

CASE  4 

CASE  5 

CASE  6 

CASE  7 

CASE  8 

CASE  9 

CASE  10 

1 

o.o 

0,129570 

0.015231 

0.006761 

0.02*991 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2 

O.IQOJUO 

0. 062*80 

0.009628 

0.004525 

0.041*40 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

J 

0.200000 

0.036112 

0.007270 

0.003624 

0.042261 

0,0 

0.0 

0.0 

0.0 

0.0 

0.0 

A 

0.400000 

0.016367 

0.  006  730 

0.004280 

0.  044040 

0*0 

0,0 

0,0 

0.0 

0.0 

0.0 

5 

0.  <500000 

0.00 7*76 

0.016049 

0.005106 

0. 033779 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

22 

i. 000000 

0.305186 

0.058976 

0.003649 

0.019742 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2J 

1. 094909 

0. 303300 

0.009163 

0.002402 

0. 008359 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

24 

1.500000 

0.001203 

0.001247 

0.000538 

0.002077 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

25 

3.000000 

0.000  202 

0.000111 

0.000026 

0.000261 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

DETA  1 LED 

LEADING  EDGE 

LOADING 

1 

0. 020000 

0.168371 

0.022600 

0.010083 

0.047722 

2 

0.04*0000 

0.127014 

0.015845 

0.007136 

0. 039566 

5 

0. 040000 

0.097  176 

0.012793 

0.  005  832 

0.038460 

4 

O.OHOOwO 

0.0’7456 

0.010931 

0.005047 

0. 039619 

5 

0.  100000 

0.062480 

0.009628 

0.00*529 

0.041840 

* CH3R0WISE  LOADING  FOR  AU  FUNDAMENTAL  CASES  • 


SECTION 

2 Y » 

0. 887500 

CHORO  ■ 

0.444444 

MING 

1 

XB 

CASE  l 

CASE  2 

CASE  3 

CASE  4 

CASE  5 

CASE  6 

C ASF  7 

CASE  8 

CASE  9 

CASE  10 

6 

0.0 

0.185800 

0.024182 

0.010750 

0.056734 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

7 

0.100000 

0.107434 

0.0I7Q9* 

0.008019 

0.  068596 

0.0 

0.0 

0.0 

O.G 

0.0 

0.0 

6 

0.200000 

U.  064708 

0.013036 

0.006439 

0.  06  7 8 74 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

. 

9 

C.  500000 

3. 0 2923t» 

0. 012062 

0.007998 

0.070135 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 1 

10 

0.900000 

0.0 12979 

0.027342 

0.011413 

0.051704 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

JET 

26 

1 .000000 

0.009876 

0.071214 

0.008014 

0.030824 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

27 

1.093999 

0.005497 

0.015626 

0.004878 

0.014111 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

I 

28 

1.50O000 

0.J01876 

0. 001991 

0.000959 

0.003273 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

J 

29 

3.000000 

0. 000287 

0. 00015? 

0.000032 

0.000369 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

detailed  leading  edge  loading 


1 

0.020000 

0.2  7368  7 

0.036243 

0.016196 

0.090728 

2 

0.043000 

0.189277 

0.  02  5879 

0.011673 

0.072112 

3 

0.060000 

0.150051 

0.021396 

0. 009764 

0.067293 

4 

0.080000 

0.125341 

0.018803 

0.0086  9* 

0.  066906 

5 

0. 1000O0 

0. 10743* 

0.017094 

0.008019 

0.068596 

I 


* CHOROMISE  LOADING  FOR  All 

FUNDAMENTAL 

CASES  • 

‘r 

SECT  1DN 

3 Y ■ 

0.687500 

CHORD  - 

0.444444 

MING 

1 

XR 

CASE  l 

CASE  2 

CASE  3 

CASE  4 

CASE  5 

CASE  6 

CASE  7 

CASE  8 

CASE  9 

CASE  10 

11 

0.0 

0.2  30761 

0.03*686 

0.01*448 

0.091838 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

12 

0.100000 

0.  1*0  8°  i1 

0.025102 

0.011127 

0.095908 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

13 

0.200000 

u. 05958* 

0.019423 

0.009298 

0. 090418 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

14 

0. 500000 

0.0*4098 

0. 01*778 

0.012470 

0.089254 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

15 

0.8000  00 

0.02514* 

0.027830 

0.032386 

0.073930 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

16 

0.90O000 

0.0 2002 5 

0.035635 

0. 085626 

0.062237 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

JET 

30 

1 .000000 

0.0  13414 

0.076381 

0.026293 

0.039217 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

31 

1 .099995 

0.UO8238 

0.019171 

0.011354 

0.019388 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

32 

1.500000 

0.00267* 

0.002923 

0.001  1 75 

0.  00*705 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

33 

3.000000 

0.000348 

0.0001*9 

0.000040 

0.00044 2 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

DETA 1LE0 

LEADING  EDGE  LOADING 

1 

0.020000 

0 . 141*09 

0.052102 

0.021837 

0.  143R41 

2 

0.040000 

0.238066 

0.037353 

0. 0 15  82 7 

0. 1 10679 

3 

0.060000 

u.  190841 

0.0310*0 

0. 01333? 

0. 099851 

4 

O.OBOOUO 

0.  lol645 

0.027*36 

0. 01 I e65 

0.  096  1 98 

5 

0.  IOOOOU 

0.  14099H 

0.025102 

0.011127 

0.  095008 

• CHC/PDMlSf  LOADING  FOR  ALL  FUNDAMENTAL  CASES  • 


SECT  I3N 

4 Y « 

0.2  75000 

CHORD  ■ 

J. 44444* 

1 

*8 

CASE  1 

CASE  2 

CASE  3 

CASE  4 

CASE  5 

CASE  6 

CASE  7 

CASE  8 

CASE  9 

C At 

17 

0.0 

0.26 1299 

0.04*340 

0.011*30 

0.121182 

0.0 

0.0 

D.D 

D.D 

0.0 

0.0 

Id 

0.  IJOOOC 

0.142285 

0. 032007 

0.007930 

0.117146 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

19 

0. 200000 

0.10520c 

0.  0?*>658 

0.005608 

0.106311 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

23 

0.  *>00000 

0.054223 

0. 022370 

0.004253 

0. 101691 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

21 

0.900000 

0.025  307 

0.036868 

0.002913 

0.071827 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

34 

l.OOOOUO 

0.117243 

0. 07*280 

0.002  095 

0.0*4616 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

35 

1.099999 

0.  ) 106b*. 

0.020564 

0.001*40 

0.022860 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

36 

1 .500000 

0. 303524 

0.00  35  78 

0. 0006*4 

0.005991 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

37 

3.000000 

u.  0(10*0  i 

0.  000220 

0. 000081 

0. 000522 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

DETA  Ilff) 

LEADING  EOGk  LOADING 

1 

0.020000 

0.337136 

0.06  6 587 

0.017101 

0.18846? 

2 

C .040 ijO 

0.27066* 

0.0*7*17 

0.01217? 

0.  142595 

1 

0.0b 00 00 

u. 2177*1 

0. 039630 

0.01002  3 

0.  126244 

4 

0.0*0000 

C. 18522" 

0. 035007 

0.00**o7 

0.  119*95 

5 

0.  100000 

0 • 1 4 2 2*  5 

0. 032007 

0. 007c  30 

0.  117146 

A5 
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V,  : *- 


«•  SPANMISE  LOAOING  FOB  FUNDAMENTAL  CASE  I • 


G.475J0C  O.'j3*»204  0. 017*53  0.054737  • 

0 . 867500  0.U5M86  0.017453  0.073639  * 

U . 08 7 *>oO  0.076014  0.017453  0.093468  • 

0.275000  0 .05  ■*  566  0.047453  0.407020  • 

--------  — • 

TCTAL  0.U7-1903  0.017453  0.096356  • 


COG 

0.0005904 
0.0009806 
0.0013267 
0 .00  456  32 


INOUCEO  OBAG 
CO 


0.0004523  0.0002930  0.0004577  1.0000000  0.0001921  0.0427911 

0.0001523  0.0006025  0.0005304  ». 0000000  0.0133976  0.0210746 

0.0004523  0.0009294  0.0005496  1.0000000  0.0101900  0.0120300 

0.0001523  0.0041917  0.0005239  1.0000000  0.0215096  0.0065474 


P ITCHING  MOMENT 


Y CMr,  CMMU  CMT  CM 

J.975D0U  -0.007645  -0.017453  0.017453  -0.007645 

0.r8T5u0  -0.013358  -0.017453  0.017453  -0.043388 

0.68750u  -0.0  l‘>  39  6 -0.017453  0.017453  -0.  019396 

0.275 OOu  -0.023761  -0.017453  0.017453  -0.023761 

TOTAL  -0.020468  -0.017453  0.017453  -0.020458 

-0.000732  -0.013090  0.013090  -0.000732 


0.0001523  0.0010010  0.0005205  0.9999996 

...  LIFT  CENTEB  ... 

• • XCP/C  XCl/C 

• • 0.222905  0.405106 

• • 0.230204  0.410010 

• • 0.255167  0.394250 

• • 0.265205  0.305100 


I APE  X)  0.259282  0.393450  IX/CRfFI 

IXMC1  0.115236  0.174007  IX/B/21 


• SPANHlSE  LOADING  FOB  FUNDAMENTAL  CASE  2 * 


0. 974 Jwo  0.ul»256  0.017463  0.030710 

0.007500  0.021353  G.017h53  0.038806 

O.6P750O  0.021*46  0.017*63  0.045009 

0.275000  0. J 3 362 7 0.01?453  0.0509E0 

TOTAL  0.020594  0.017403  0.047047 


COG 

0.0 

0.0 

0.0 

0.0 

CDMU 

0.0001523 

0.0001523 

0.0001523 

0.0001523 

CS 

0.0000040 

0.0000102 

0.0000210 

0.0000343 

CO 

0.0001483 
0.0001421 
0.3001313 
0.000 1160 

CMU 

1.0000000 

1 .9  000000 
1.0000000 
1.0000000 

GAMMA 

0.0044180 

0.0068049 

0.0009210 

0.0103034 

ALFIN 
0.0254144 
0.0122251 
0.0057449 
0.002  2070 

0 .0 

0.0001523 

0.0000261 

0.0001262 

0.9999996 

0.0001164 

PITCHING  MQMFNT 


Y CMC  CMMU  CMT 

0.  r>7500o  -0.007630  -0.017453  0.0 

0 . P07*  00  -0.012071  -0.017*43  0.0 

O.«076jo  -0.01414*  -0.017*53  0.0 

0./75J0O  -J.J170-H  -0.017*53  0.0 

TllTAl  -0.014**0  -0.017453  0.0 

-0.00*0*.!  -O.O13040  0.0 


CM 

-0.  026203 
-0.  029524 
-0.032547 
-0.  034498 


-0.032843  I A PE  X I 
-0.021131  IXMCI 


...  LIFT  CENTEB  ... 

XCP/C  XCL/C 
0.590655  0.823296 

0.5  6 5 3 05  0.7  60  0 1 2 

0.532393  0.710203 

0.500309  0.676695 

0.521710  0.699142  IX/CBEFI 

0.231071  0.310730  IX/0/21 


• SPANMISE  LGAIING  FOR  FUNDAMENTAL  CASE  3 • 


I IF? 

y cu.  r.iMvi 

O.C75JOO  0.00*011  0.0 

0.8*75  JU  0.JO-:'23O  0.0 


CL 

0.004811 

0.009230 


0 .617*.  00  0.02*410  0.011453  0.041964  * 
0.275000  0.005419  0.0  0.005419  • 


0.011115  0.00*800  0.315915  • 


COG 

0.0 

0.0 

0.0000977 

0.0 


INDUCED  DRAG 

CS  CD 

0.0000000  -0.0000000 


CMU  GAMMA  ALFIN 

1.0000000  0.0013792  0.0050237 


0.0  0.0000020  -0.0000020  1.0000000  0.0026257  0.0014037 

0.0001523  0.0000036  0.3002463  1.0000000  0.0066333  0.0100795 

0.0  0.0000023  -0.0000023  1.0000000  0.0014751  0.0095400 


0.0000419  0.0000025  0.3000662  0.9999996 


v C MG  CMMU 

(I. *: 75 300  -0.03224P  0.0 

0.f876Jj  -0.00*600  0.0 

0. *87500  -0.014220  -0.017*53 
0.275300  -0.001887  0.0 


P ITCHING  MOMENT 
CMMU  CMT 

0.0  0.0 

0.0  0.0 


CMT  C M 

0.0  -0.00224S 

0.0  -0.004600 

0.001745  -0.031936 
0.0  -0.001  867 


XCP/C 

0.467135 


XCL/C 

0.467105 


TOTAL  -0.006189  -0.00*800  0.0004EO  -0.010509 

-0.003410  -0.003600  0.000*80  -0.006530 


0.499208  0.499200 

0.662069  0.802619 

0.343260  0.348260 

0.556798  0.690460  (X/CREFI 

0.247466  0.306071  (X/B/2) 


SPAMWISE  LCADING  F Ok  FUNDAMENT  A|  CASE  4 * 


Y CLG  Cl  Mu  CL  * 

0.975 JOo  0.040381  0.026180  0.066561  • 

0.067500  0.065122  0.02*180  0.091302  • 

0.6875 00  u . u8?»  40*  0.02j180  0.112785  • 

0.27500U  0.100)1*  0.026 180  0.12709*  * 


INDUCED  ORAG 
CD 


0.0003055  0.0003427  0.0000147  0.0006335  1.0000000  0.0101972  0.0550330 

0.000*605  0.0003427  0.0000562  0.3007550  1.0000000  0.0164325  0.0277257 

0.0005723  0.0003427  0.0001472  0.00076?e  1.0000000  0.0219236  0.0144727 

0.0006315  0.0003427  0.0002580  0.0007162  1.0000000  0.0256589  0.0069044 


3.009*78  0.026  IDO  0.1  15658  • 0.0005785  0.0003427  0.0001901  0.0007311  0.9999996 


PITCHING  MOMENT 


...  LIFT  CENTER  ... 


V CMC  CMMU 

CMT 

CM 

* 

XCP/C 

XCL/C 

0.  )75JOO  -0.019607  -0.026100 

0.0 

-0.044787 

• 

0.4607  76 

0.672465 

0.887500  -0.029265  -0.026100 

0.0 

-0.055445 

* 

0.449387 

0.607270 

3.68750O  -0.037602  -0,026100 

0.0 

-0.063702 

* 

0.434162 

0.565522 

0.27500O  -0.042400  -0.026180 

0.0 

-0.068670 

• 

0.421055 

0.540310 

TOTAL  -0.038299  -0.026180 

0.0 

-0.064*79 

1 APE  X) 

0.42P023 

0.557493 

IX/CREF) 

-0.01592)  -0.019635 

0.0 

-0.035564 

IXMCI 

0.190232 

0.247775 

IX/ 8/21 

I I FT  COEFFICIENT  DERIVATIVE  DUE  TO  PITCHING  ABOUT  XCG.  CLO  ■ 0.009470 

PITCHING  MOMENT  COEFFICIENT  DERIVATIVE  AB'JUT  ORIGIN  OUE  TO  PITCHING  ABOUT  XCG.  CMQ  • -0.030299 
PITCHING  MOMENT  CO  FF  F DEPlVATIVE  AROUT  XMC  OUE  TO  PITCHING  ABOUT  XCG.  CMQMC  ■ -0.015929 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
FRUJM  COtY  EUKNISHED  TO  DDC " 


* « « 


* TOTAL  AERUQYKAKIC  COEFFICIENTS  * 


CASE  1 

CASE  2 

CASE  3 

CASE  4 

CASE  5 

CASE  6 

CASE  7 

CASE  0 

CASE  9 

CAS 

CC  LG 

O.O704O2E 

0.02959*1 

0.0111152 

0.0896784 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CCLJ 

0 *0  1 7 AS  9 3 

0.017*533 

0.00*7997 

0.0261799 

0.0 

0.0 

9.0 

0.0 

0.0 

0.0 

••  CC  L 

0.056  3561 

0.0*  70*7* 

0.0 1591*9 

0.1156553 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CCC6 

0.331 3771 

0.0 

0.0000269 

0.0005785 

0.0 

0.0 

0.3 

0.0 

0.0 

0.0 

CCOJ 

0.300  152  3 

0.0001523 

0.0000419 

0.0003*27 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CCS 

0.00  10013 

0. 3000  2o 1 

0.0000025 

0. 0001 901 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

ccn 

0,0009295 

0.0001262 

0.0000662 

0.0007311 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

••  coir/ 

0,00010 33 

0.300116* 

3.0000785 

0.0007183 

0.0 

0.0 

0.0 

0.0 

3.0 

0.0 

••  CC  J 

0.99<*V99F 

0.9999996 

0.9999996 

0.9999996 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CCMG 

-0.020*581 

-0.315*396 

-0.0061859 

-0.03F  2480 

0.0 

0.  0 

0.0 

0.0 

0.0 

0.0 

CC  MJ 

-0.017*533 

-0.017*533 

-0.0047997 

-0.02 t 1799 

0.0 

0.0 

0.9 

0.0 

0.0 

0.0 

CC*T 

0.017*533 

0.0 

0.0004800 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CC* 

-0.020*581 

-0.0328923 

-0.0  1050  86 

-0.06**787 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CXCP 

0.2*92817 

0.5217098 

0.5*67951 

0.4280229 

0.0 

0.0 

0.0 

0.0 

3.0 

0.0 

CXCL 

0.393*50* 

0.6991*21 

0.6904603 

0.5574928 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CXCPB 

0.1152363 

0.2318709 

0.2*74650 

0.1902324 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CxCLB 

0.1758668 

0.310729E 

0.3068712 

0.247  7745 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CCHGMC 

-0.000  732* 

-0.0080*10 

-0.0034101 

-0.0159292 

o.o 

0.0 

0.0 

0.0 

0.0 

0.0 

CCMjHC 

-0.0130900 

-0, 01 30*00 

-0.0035997 

-0.0196349 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CLNTMC 

0.0130930 

G.O 

0.0004800 

0.0 

0.0 

0.0 

3.0 

0.0 

0.0 

0.0 

•*  CCMMC 

-0.000732* 

-0.0211 310 

-0.0065299 

-0.035564 2 

0.0 

0.0 

3.0 

0.0 

0.0 

0.0 

CLLG 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CILJ 

0.3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

• CLl 

0.0 

0.0 

0.0 

0,0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CNJ 

0 .0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

• CNIMC 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

• CC  V 

0.0 

0.3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CBGF 

0.0358229 

0.01346*0 

0.0067121 

0.0409301 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CBGL 

0.0358229 

0.013*6*0 

0.006  7121 

0.04  03301 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CBJP 

0 .008726* 

0.0057266 

0.0032998 

0.0130900 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CBJL 

0.0087266 

C. 0087266 

0.0032998 

0.0133900 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

ce* 

0.0**5495 

0.3221906 

0.01001 19 

0.0539200 

0.0 

0.0 

3.0 

0.0 

0.0 

0.0 

CBL 

0.0** 5*95 

0.0221906 

0.0100119 

0.0539200 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CPMBF 

0.462  3*23 

0.47166*7 

0,  *290900 

0.4662000 

0.0 

0.0 

0.0 

0.0 

9.0 

0.0 

CPHBl 

0. *623*23 

O.*7lo6*7 

0.6290900 

0.4662009 

3.0 

0.0 

0.0 

0.0 

0.0 

0.0 

THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
FROM  COP1  FURNISHED  TO  DDC . 


A7 


• CHOIOMlSE  l PAD  INC  FOR  COMPOSITE  CASE  1 • 


FUNDAMENTAL  CASE  FACTORS 

A<*»  Al  3 • A (A  I A ( 51  AI6)  AI7»  AI8)  AI9)  AIIOI 

0.0  10.000000  10.000000  0.0  0.0  0.0  0.0  0.0  0.0  0.0 

**•  A.OTE  EACH  LEADING  EDGE  CP  VALUE  IS  THE  AVERAGE  VALUE  OF  THE  SINGULAR  DISTRIBUTION 

OU  NIT  PLOT  THESE  LOAOlNG  POINTS  DIRECTLY 

SECTION  1 V ■ 0.975000  CHORD  ■ 0. ****** 

MING 

KB  0.0  0.100000  0.200000  0.500000  0.900000 

CP|A>0)  0.219527  0.1*156*  0.1099*1  0.110103  0.211553 

CPIA-ll  0.129570  0.062*80  0.036112  0.016367  0.007*76 

JET 

XB  1.000000  1.099999  1.500000  3.000000 

CP  I AO  I 0.6262**  0.  115851  0.018*50  0.001371 

CPIA-ll  0.005186  0.003300  0.001203  0.000202 


SECTION  2 

MING 

XR  0.0  0. 100000  0.200000 

CPIA-OI  0.3*9325  0.251132  0.1953*5 

CPIA-ll  0.18580  0 0.107*1*  0.06  *7  08 

JET 

xe  1.000000  1 .09  .999  1.500000 

CPIA-OI  0.752272  0.205033  0.029*91 

CPI  A- 1 ) U. 008878  0.005*57  0.001876 

SECTION  3 Y - 0.687500  CHOPD  ■ 0.****** 

0.0  0.100000  0.200000  0.500000  0.000000  0.900000 

0. *913*0  J. 362289  0.7P7105  0.312*86  0.602161  1.212606 

u. 230761  0.1*0898  0.08958*  0.0**098  0.0251*7  0.020025 

1.000000  1.099999  1.500000  3.000000 

1.0267*3  0.30525 2 0.0*0979  0.002192 

0.017*1*  0.008238  0.00267*  C.0003*3 


SECT I3N  * v • 0.275000  CHORD  ■ 0.4***** 

MING 

XH  0.0  0.100000  0.200000  C. 500000  0.900000 

CPIA-OI  0.557696  0.399370  0.30*666  0.26623*  0.397805 

CPIA-ll  0.2612°9  0.162285  0.105206  0.05*223  0.025307 

JET 

xe  1.300000  1.099Q99  1.500000  3.000000 

CPIA-OI  0.d0J7**  0.2200*1  0.0*2222  0.003008 

CPIA-ll  0.0172*3  0.010685  0.00352*  0.000*03 


* COMPOSITE  CASE  1 • 


FUNDAMENTAL 

CASE  FACTORS 

All) 

A<  21 

Al  31 

Al*  1 

A ( 51 

A 1 6 ) A ( 7 1 

Al  8) 

A 1 9 ) 

Al  10) 

3.0 

lu.o.ioooo 

10.  JOOoOU 

0.0 

0.0 

0.0  0.0 

0.3 

0.0 

0.0 

SECTION 

Y 

C LOO 

CLMUO 

CLO 

* 

CMGO  CMMUO 

CMTO 

CNO 

• 

XCPO/C 

x:lo/c 

CL‘»A 

CLMUA 

CLA 

• 

CMGA  CMMUA 

CMT  A 

CNA 

• 

XCPA/C 

XCL  A/C 

1 

0.97 5000 

0. 1M067* 

0. 17*533 

0.355206 

• 

-0.100776  -0.17*533 

0.0 

-0.275309 

• 

0.557779 

0.775067 

0.3  1*2  8** 

0.01 7*53 

0.051737 

• 

-0.0076*5  -0.017*53 

0.017*53 

-0.0076*5 

• 

0.222985 

0. *85106 

2 

0.887500 

0.3)5825 

0. 17*533 

0. *80358 

• 

-0.16676*  -0.17*533 

0.0 

-0.3*1317 

• 

0.5*5357 

0.7105*7 

0.3*61  86 

0. 0 1 T*5  3 

0.073639 

• 

-0.013388  -0.017*53 

0.017*53 

-0.013388 

• 

0.23828* 

0. *18818 

3 

0.687500 

0.5.!9560 

0.  3*9066 

0.87*626 

• 

-0.313718  -0.3*9066 

0.017*53 

-0.6*5331 

6 

0.592*13 

0.75*3*2 

0.0  '60  1* 

0.017*53 

0. 093*68 

• 

-0.019396  -0.017*53 

0.317*53 

-0.019396 

• 

0.255167 

0.39*250 

* 

0. 275000 

0.369*60 

0.17*533 

0. 563993 

• 

-0.189319  -0.17*533 

0.0 

-0.363852 

• 

0.  *86  1 08 

0.6*5136 

0.0«  >566 

0.017*53 

0.  107020 

• 

-0.023761  -0.017*53 

0.017*53 

-0.023761 

• 

0.265285 

0.385106 

TOTAL 

0.4  » TO 93 

0.222529 

0.629623 

* 

-0.216285  -0.222529 

0.00*800 

-0.*3*0 1* 

1 APE 

X) 

0.531290 

0.6969*8 

0.0' 3903 

0.017*53 

0.096356 

• 

-0.020*58  -0.017*53 

0.017*53 

-0.020*58 

I APE 

XI 

0.259282 

0.393*50 

m 

-0.11*511  -0.166897 

0.00*800 

-0.276609 

IXMC  1 

0.236129 

0.30975* 

• 

-0.000732  -0.013090 

0.01 3090 

-0.000732 

IXMC  ) 

0.115236 

0.17*867 

SI  C T If  N 

1 

CDGO 

C09U3 

CSO 

...  INDUCED 
COO 

DRAG 

G A MM AO 

AL  F "l  NO 

CTO 

CNU 

CDGX 

CDMUX 

CSX 

CDX 

COGA  2 

CDMUA2 

CSA2 

COA  2 

GAMMA2 

ALE  INA 

1 

0.97?3JJ 

0.0 

0.01*2  109 

0.00084*2 

0.0143867 

0.0579720 

0.30*3813 

0.9856133 

1.0000000 

0.3031*33 

0.0030*62 

0.0009947 

0. 00 520* f 

0.000*984 

0.0001523 

0 . 00029  30 

0. 000*577 

0.008 1921 

0.0427911 

2 

0 .ReVOO 

0.0 

0 • Cl  52  309 

0.0021298 

0.0131011 

0.095 1066 

0.1370882 

0.98  6 8 9 9 0 

1.0000000 

0.0051)77 

0.00  30*62 

0.00226  56 

0.0061182 

0.00)9.106 

0.0001523 

0.0006025 

0.000530* 

0.0133976 

0.02 187*6 

3 

0 .6*  75 JO 

0.0195*20 

0.06 0°2  35 

0.00*21  35 

0.0762519 

0.1555502 

0.1662**3 

0.9237*81 

1.0000000 

0 .Uo9*  )4«. 

0.0060923 

0.0039578 

0.0116690 

0 .031)267 

0.0001*23 

0. 000920* 

0.0005496 

0.0181988 

0.0120300 

* 

0.2  7‘OJO 

0.3 

0.01  52)0*9 

0.005*28* 

0.00  9802* 

0.11778*9 

0.117*700 

0.9901975 

1.0000000 

0.006  *9*4 

0.00  30*6 2 

0. 0050868 

0.00*7568 

0.001*6  )2 

0.0001523 

0.001 1917 

0.0005239 

0.0215896 

0.0065*74 

TOTAL 

3 .30*  )T«.l 

0.0277963 

0.0044528 

0.0287176 

0.3451221 

0.9712820 

0.9999996 

0. J3716  *4 

0. 00)*  8 3° 

0.3042190 

0.0068502 

0.3101239 

J.0tl|)77| 

0.0001523 

0.0010010 

0.0005285 

0.3005033 

MIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
A8  IRON  COlrY  EURtilSttFD  TO  DDG  . — ' 


* T3TAL  AEAOOYNANIC  COCMIC1FNTS  • 


ALPHA  *0 

CCLG 

0. *070*3 

CCLJ 

0.22252* 

*•  CfL 

0.629623 

CCOG 

0.0053  7*1 

CCOJ 

0.0277963 

CCS 

0.00**526 

CCD 

0.0287176 

••  COIW 

0.0*51221 

••  CCJ 

0.9**9996 

CCHG 

-0.216285 

CCMJ 

-0.222529 

CC**T 

0.00*800 

CC* 

-0.*3401* 

:«cp 

0.5312*0 

CMCL 

0.6969*8 

CKCPb 

0.23612* 

CHCLB 

0.30*75* 

CC«G**C 

-0.11*511 

CCMJHC 

-0.166897 

CCMTMC 

0 .00*800 

••  CCMHC 

-0.27660* 

CLLG 

0.0 

CLLJ 

0.0 

* CLL 

0.0 

l.NJ 

0.0 

• CNIHC 

0.0 

• CCV 

0.0 

CBGR 

0.201761 

CBGL 

0.201761 

CBJR 

0.1202  6* 

C8JL 

0.12026* 

CBR 

0.122025 

C8L 

0.322025 

CPMhR 

0.511*57 

CPHBL 

0.511*57 

ALPHA  ALPHA**? 

0.078*01 
0.017*53 
0.0*635* 

0.007185*  0.0013771 

0.001881*  0.00015?! 

0.00*21*0  0.0010010 

0.0068502  0.0005285 

0.010121*  0.0005011 

- 0. 020*58 
-0.017*53 
0.017*53 
-0.  020*58 
0.25*28? 

0.393*50 
0.115236 
0.17*867 
-0.000712 
-0.0130*0 
0.013090 
-0.  000732 
0.0 
0.0 
0.0 

0.0  0.0 

0.0  0.0 

0.035823 
0.035823 
0.008727 
0.008727 
0.0**5** 

0.  0**5*9 
0. *623*2 
0. *623*2 
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• SECOND  RUN  FOR  STABILITY  DERI  VAT | VE  CASE  * 


• STABILITY  DERIVATIVE  DATA  FOR  FUNDAMENTAL  CASE  I • 


POLLING  MOMENT  COEFFICIENT  DERIVATIVE  DUE  TO  YAWING*  CLLR  - 0.0001160 

YAWING  MOMfur  COEFFICIENT  ABOUT  XMC  OUE  TO  YAWING  ABOUT  XCG,  CNIRI  MAY  BE  CALCULATED  AS  FOLLOWS 

CNIRI  * CNR*R  ♦ CNR2*R«*2 

WHERE  CNR  • 0.00000 1049 

CNR  2 • 0.0 

SIDE  FORCE  COEFFICIENT  DUE  TO  VAwiNG.  CVIFI  MAY  BE  CALCULATED  AS  FOLLOWS 

CYIRI  * CVR*R  ♦ CYR2*R**2 

WHERE  C VR  ■ 0.0 

CVR2  • 0.0 


• STABILITY  DERIVATIVE  DATA  FOR  FUNDAMENTAL  CASE  2 • 


MILLING  MOMENT  COEFFICIENT  DERIVATIVE  OUE  TO  YAWING.  CLLR  • 0.0000436 

YAWING  MOMENT  COEFFICIENT  ABOUT  XMC  DUE  TO  YAWING  ABOUT  XCG,  CNIRI  MAY  BE  CALCULATED  AS  FOLLOWS 

CNIRI  * CNR*R  ♦ CNR2*R**2 

WHERE  CNR  « 0.000000034 

CNR 2 - 0.0 

SIDE  FORCE.  COEFFICIENT  DUE  TO  YAWING.  CYIRI  MAY  BE  CALCULATED  AS  FOLLOWS 

CYIRI  - C YR* R ♦ C YR2*R**2 

WHERE  CYR  - 0.0 

CYR2  « 0.0 


• STABILITY  DERIVATIVE  DATA  FOR  FUNDAMENTAL  CASE  3 • 


MILLING  MOMFNT  COEFFICIENT  DERIVATIVE  OUE  TO  YAWING.  CLLR  * 0.0000344 

YAWING  MOMfc^T  COEFFICIENT  ABOUT  XMC  DUE  TO  YAWING  ABOUT  XCG,  CNIRI  MAY  BE  CALCULATED  AS  FOLLOWS 

CNIRI  « CNR4R  * C NR2* R**2 

WHERE  CNR  ■ -0.000000100 

CNR2  ■ 0.0 

SIDE  FORCE  COEFFICIENT  DUE  TO  YAWING,  CY(R|  MAY  BE  CALCULATED  AS  FOLLOWS 

CYIRI  ■ CYR*P  ♦ CYR2*R**2 

WHfRE  C YR  » 0.0 

CYP2  • 0.0 


♦ STABILITY  DERIVATIVE  DATA  FOR  FUNDAMENTAL  CASE  4 • 


POLLING  MOMENT  COEFF  DERIVATIVE  DUE  TC  ROLLING,  CLLP  ■ -0.0066944 

YAWING  MJMENT  COEFFICIENT  AHDUT  XMC  DUE  TO  POLLING.  CN(P|  MAY  BE  CALCULATED  AS  FOLLOWS 

CNIP1  ■ CNP24P**2 
WHERE  CNP2  ■ 0.0 

SIDE  rciRCF  COEFFICIENT  DUE  TO  ROLLING,  CVIPI  MAY  BE  CALCULATED  AS  FOLLOWS 
CYIPI  ■ CYP2 *P**2 


WHERE  CYP2  ■ 0.0 


jMM  CQtl  FUB»ISHH>  TO  — 


* t *"•"*? - i 
. - ..i  ’ :*■  * . . 


- STABILITY  DERIVATIVE  DATA  FLR  COMPOSITE  CASE  1 • 


FUNDAMENTAL  CASE  FACTO** 

A||»  At  21  All)  A 1 4 ) A 1 51  AI6I  At  7 » Alt)  A(*l  A<  101 

0.0  <0.0)0000  10.000000  0.0  0.0  0.0  0.0  0.0  0.0  0.0 


LIFT  COLFriCIENT  DERIVATIVE  OUE  TO  PITCHING  ABOUT  ICG  . CLO  • 0.0094TB 

PITCHING  MOMENT  COEFFICIENT  DERIVATIVE  ABOUT  ORIGIN  OUE  TO  PITCHING  ABOUT  ICG.  CNQ  • -0.0JB299 
PITCHING  MOMENT  COFFF  DERIVATIVE  ABOUT  INC  OUE  TO  PITCHING  ABOUT  ICG.  CMQMC  - >0.015525 


ROLI  1NG  MOMENT  COEFF  DERIVATIVE  OUE  TO  ROLLING.  CUP  ■ -0.0046549 


Y AMI  KG  MOMENT  COEFFICIENT  ABOUT  KMC  OUE  TO  ROLLING.  CNCPI  MAY  BE  CALCULATED  AS  FOLLOWS 
CNIP)  • CNP9P  ♦ CNP2*P**? 

WHERF  CNP  • CNPO  ♦ CNPA*ALPHA 
CNPO  ■ -0.000.T5Jl 

CNPA  • -0.0000601 

CNP?  - 0.0 


SIDE  FORCE  COEFFICIENT  OUf  TO  ROLLING.  CVCPI  MAY  BE  CALCUIATEC  AS  FOLLOWS 
C V I P ) ■ CVPPP  ♦ CVP2*P**2 
WHERE  CVP  • CYPO  ♦ C YPA*Al PHA 

CVPO  ■ 0.0 

CYPA  - 0.0 

CYP?  - 0.0 


RL'LL  INC  MOMENT  COEFF  DERIVATIVE  DUf  TO  YAWING  ABOUT  XCG . CLLR  MAY  RE  CALCULATED  AS  FOLLOWS 

CLLR  > CLLRO  A CLLR A9ALPHA 
WHERE  CLLRO  - 0.000700$ 

CLLRA  - 0.0001I6B 


APPENDIX  B - ERROR  MESSAGES 
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In  order  to  avoid  wasting  computing  time,  some  types  of  errors  or 
inconsistencies  in  input  data  can  be  checked  by  the  program  before  beginning 
computation.  Sometimes  such  errors  make  the  program  unable  to  continue  execu- 
tion, and  sometimes  the  program  can  recover  by  ignoring  the  data  or  substi- 
tuting correct  data.  In  any  case,  an  error  message  will  be  printed,  enabling 
the  user  to  identify  the  type  of  error  or  warning  him  that  the  input  data  may 
have  been  altered.  Below  is  a list  of  the  messages  which  may  be  printed  and 
the  subroutines  in  which  they  occur.  At  the  end  of  each  run  a status  message 
is  always  printed  indicating  whether  the  program  has  reached  a normal  or 
abnormal  conclusion. 


Message 

-r-*  *«.*«!*  » -I 

* T,«r  PROGRAM  HAS  KT.  ACMEI)  TERMINATION  * 

• T Ht  PPOoPA*  HA  S RFACmH)  ABNORMAL  TERMINATION  • 

fSY***  H '•  j IWlCAtCO  AA  ANT  J-SY/'f  TVJC  CASE*  HOwE VI R » IT  WILL  BE  TREATED  AS  SYNEIRIC. 

fMfc  Hi  TV  PC  FI  Ao  MAS  IHF  VALUE  OF** 

ONLY  THE  VALUES  l OR  2 A«F  ACCEPTABLE 

THIS  CASE  MAS  TEEN  TERMINATED 

KROWS  ■»•* 

MJfMFR  ME  WING  Row  TYPES  ■♦** 

*♦*  WING  FLFMrNTS  PRESCRIBED  FOR  ROW  TVPF*** 

NUMBER  TE  JCT  ROW  TYPES 

*^«  JCT  t-U."ENTS  PRESCRIBED  FOR  HOW  TYPC*** 

3 ROW  Cnt,nn>)ltV  PULL  FAILURE 

A'l  1'ic  '.‘ISIS  l-NCY  N-xS  F UNO  IN  TMF  SECTIONAL  LtAOlNG  ANO  TRAILING  EDGE  INPUT 

**L  F AS  * LHCl.K  YOllR  SFClION  LCiCATION  < V I INPUT 
***♦  IS  Tf.r  MANY  ELEMENTS 

A /PRO  VALUE  OF  C M||  HAS  BEEN  INPUT.  THIS  CMU  CASE  HAS  BEEN  IGNORED 


Subroutine 


STAGE  I 
I .VP  T S 

INPTS 


IVP'JT  J 
INPUT  J 


INPIITJ 


XLETRl 

BOAS 


MUXS 

noxj 


MJN TA.NENTAL  GFUME  TR  1C  CASE*** 

AM  IMCOMsf  Sir  UCY  HAS  'J*FN  t-C  tJNI  IN  THE  HINGE  INPUT  DATA  FOR  WING  ROW***#  ROW  TYPE*** 


vu  INC.  P*acT  CMMP'  SJTf  (.4 SF  INMUI  VALUE  HAS  BEEN  FOUNO.  IT  WILL  RE  IGNORED. 
HHl  IMAM  ’>  CO- TP'  IS  I T F CASTS  HAVE  REIN  INPUT.  SUBSEQUENT  INPUTS  WILL  BE  IGNORED. 


INCOMP 

INCOMP 


v.  IMP  ' Flit  MAS  RIFT.  P *■  Ai)  TOPING  COMPOSITE  CASE  INPUT 


'AT*'!*  )UF  S OUT  MAVfc  ENOUGH  COPT  TO  WORK 
TU|>  CASE  HAS  MtCTJ  lERMINAftO 


INCOMP 

STG2S 


B1 
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